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Abstract 



Deuteron Compton scattering below pion photoproduction threshold is con- 
sidered in the framework of the nonrelativistic diagrammatic approach with 
the Bonn OBE potential. A complete gauge-invariant set of diagrams is 
taken into account which includes resonance diagrams without and with NN- 
rescattering and diagrams with one- and two-body seagulls. The seagull oper- 
ators are analyzed in detail, and their relations with free- and bound-nucleon 
polarizabilities are discussed. It is found that both dipole and higher-order 
polarizabilities of the nucleon are required for a quantitative description of 
recent experimental data. An estimate of the isospin-averaged dipole electro- 
magnetic polarizabilities of the nucleon and the polarizabilities of the neutron 
is obtained from the data. 
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I. INTRODUCTION 



Elastic photon, or Compton, scattering is a powerful tool for probing the structure of 
hadrons and nuclei. A deformation of the system's ground state caused by an incoming 
electromagnetic wave and encoded into electromagnetic polarizabilities of the system con- 
tributes to radiation of outgoing photons and thus shows itself in such observables as the 
differential cross section of Compton scattering. A particular example is forward Compton 
scattering. The corresponding spin-averaged amplitude at sufficiently low energies u has the 
formQ 

T(u) = e • e'*( - ^ ± 4n(a + ^ + ...). (1.1) 

Here e and e' are polarizations of the initial and final photons, Ze and M t are the elec- 
tric charge and the mass of the target, and a and (3 are the electric and magnetic dipole 
polarizabilities of the target. Many efforts have been spent to measure the polarizabilities 
of the nucleon, ajv and (3n (as well as polarizabilities of other hadrons and nuclei), and to 
understand them theoretically. For a review and further references see [1] |6| . 

The polarizabilities of the proton have been successfully found in a series of experiments 
on 7p scattering |7Hl3| which ultimately yielded quite an accurate result, 

a p = 12.1 ±0.8 ±0.5, ^ = 2.1^0.8^0.5 (1.2) 



(in the units of 10 4 fm 3 used for the dipole polarizabilities throughout the paper). The 
values (|1.2j) quoted here |12| have been extracted from data of a few experiments of 90's 



performed at energies below pion photoproduction threshold under the theoretical constraint 



of the Baldin sum rule [1 1.15 



a + P 



du 



;i.3) 



where cr to t(^) is the total photoabsorption cross section. Recently, the result 
confirmed by a more comprehensive analysis of a larger data base [|ll| . 



^2l) has been 



Meantime, studies of polarizabilities of the neutron which began even earlier than those 
for the proton (see a book |17j which summarizes a long history of these studies) achieved 
a knowledge of a n and (3 n far less satisfactory. Most of the experiments performed for 
measuring the polarizabilities of the neutron had deals with neutron transmission in the 
substance. The long-range polarization interaction 



V pol (r) = --47r[aE 2 (r) + (3H 2 (r)} 



;i.4) 



1 The factor of 4-7T in Eqs. ( |l.ip and ( |1.4| ) below stands because we use Heaviside's units for the 
electric charges and electromagnetic fields (e.g., e 2 = 4-7r/137) but, for historical reasons, Gaussian 
units for the polarizabilities themselves. 
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of the neutrons with the electromagnetic (actually, electric) field near the edge of nuclei in 
the substance creates a small but detectable contribution to the total cross section due to its 
anomalous energy dependence oc \fE [r7-31| . The best results for the electric polarizability 

MM 



of the neutron, 



a,. 



obtained from these studies read 



a,. 



12.6 ±1.5 ±2.0 [20], 



0.6 ±5 P] 



;i.s) 



where a small relativistic correction = ±0.62 || missing in a fully nonrelativistic formalism 
used in the original works [T^J2yj2T|] is included.^ Since the two values in Eq. ( |1.5| ) seem to 
contradict to each other, the current situation with knowing the polarizability a n is hardly 



satisfactory. Moreover, there is an argument |23| that the systematic (in fact, theoretical) 
uncertainty, which is a very delicate problem for those experiments, may be strongly under- 
estimated in Ref. [pOfl . Furthermore, the neutron transmission experiments do not constrain 
the magnetic polarizability of the neutron at all, although /3 n can be theoretically derived 
from a n by using the Baldin sum rule ( |1.3|) . For all these reasons there is a need for searching 
for other ways of finding the polarizabilities of the neutron, e.g., by using real photons. 

There are several reasons why experimental studies of neutron Compton scattering and 
a further extraction of the neutron polarizabilities are more difficult than those for the 
proton. First, because of the absence of dense free-neutron targets, actual measurements of 
7n scattering are forced to have a deal with neutrons bound in nuclei and hence to take into 
account effects of the nuclear environment. Second, due to vanishing the neutron Thomson 
scattering amplitude (viz. the amplitude of photon scattering off the electric charge of the 
neutron which is zero), the contribution of polarizabilities of the neutron to the differential 
cross section at low energies (< 100 MeV) turns out to be rather small. It is of order 
0(u) 4 ) in the low-energy expansion over the photon energy uo vs. order 0(u 2 ) in the proton 
case. Third, the 0(u A ) contribution of the neutron polarizabilities is accompanied with 
other terms 0(uj a ) which come from the spin-dependent part of the scattering amplitude; 
these additional terms are determined by the so-called spin polarizabilities and they cannot 
be isolated in a model-independent way P^Po] . Therefore, a use of further assumptions, 

for 



like those constituting the dispersion theory of Compton scattering [p6} - p8||2^ , |29| j3( 
evaluating the unknown pieces becomes unavoidable. All that introduces larger theoretical 
uncertainties to the obtained polarizabilities which are at least ±2 even without the nuclear- 
environment corrections. 

The first attempt to measure low-energy scattering and to extract polarizabilities of 
the neutron has been done by the Gottingen-Mainz group PJJ which followed an earlier 



2 Recently, the relation between the polarizability a n standing in the 'yn scattering amplitude and 
the so-called static polarizability a n determined in the neutron transmission experiments was re- 



derived in Ref. [22|. In essence, the main conclusion of that analysis agrees with our own findings 



This agreement, however, might be difficult to see from the paper [22], because its authors 
erroneously claim that there is a difference between a, which is understood as a parameter standing 
in an effective relativistic Hamiltonian, and the polarizability a s (using their notation) determined 
through Compton scattering. 



3 



theoretical suggestion [R2] to exploit the reaction 7c? — > 'ynp in the quasi-free kinematics. 



The result of this very first experiment, 

a n = 10.7+^, (1.6) 

is not yet so accurate. However, with the use of a wider range of photon energies (up to 
200-250 MeV), further improvements are quite feasible [24,32,33]- A high accuracy of the 



underlying dispersion calculations of scattering |23]J2IJ is crucial for finding the polar- 
izabilities from Compton scattering data taken at "high" energies, viz. those above pion 
photoproduction threshold. Therefore, a determination of a n and j3 n from such data also 
assumes a careful check and normalization of pion photoproduction off the neutron which 
is used in the dispersion calculations as a crucial input. Fortunately, such a check can be 
done in parallel with Compton scattering measurements, because the 7(71, 7r) reaction can 
be learned from j(d, ir) in the quasi-free kinematics as well ||34|| . 

In the present work we analyze another possibility for probing the polarizabilities of 
the neutron which requires no strong assumptions on the "high-energy" behavior of 'jN 
interactions. This possibility mentioned already in Ref. fLM is elastic jd scattering below 



pion photoproduction threshold. The presence of the proton next to the neutron and the 
coherence of the proton and neutron contributions makes two advantages. First, the 0(uj 2 ) 
contribution of the neutron polarizabilities to the scattering amplitude can interfere with the 
(9(1) contribution from proton Thomson scattering, so that a sensitivity of the differential 
cross section with respect to the polarizabilities is enhanced. Second, the largest contribution 
to the spin polarizabilities of nucleons which comes from the t-channel 7r°-exchange does 
not contribute to 7c? scattering at all (due to isospin), so that 0(u 4 ) corrections, which 
are not small for individual nucleons, especially for the neutron, are more suppressed in 
the considered case. Nevertheless, various binding corrections, including meson-exchange 
currents (MEC) and meson-exchange seagulls (MES), are rather important and have to be 
introduced and carefully evaluated. Their analysis is the central subject of the present paper. 

Theoretical studies of deuteron Compton scattering have been started by Bethe and 



Peierls [3^] who considered this process in the dipole El approximation. After then a 
number of calculations has been performed in 50 's and 60 's, mostly based on the impulse 
approximation A higher level of art, with an explicit consideration of MEC and of their 
influence on the so-called resonance and seagull amplitudes of Compton scattering was in- 
troduced by Weyrauch and Arenhovel [B^j . They directly calculated the seagull contribution 
from the pion exchange and developed an approximate scheme based on dispersion relations 
in the long wave-length limit to find the resonance amplitude of Compton scattering from a 
theoretically known deuteron photodisintegration amplitude. Later on, direct calculations 
of the resonance amplitude free from the approximations of the oversimplified dispersion 



relations have been performed using a simple separable NN potential ||38|| . More recently, 



this consideration was further improved |39| by using realistic NN potentials for evaluation 
of NN rescattering in the intermediate state and by taking into account leading relativistic 
corrections and MEC beyond the Siegert approximation. A similar (but technically differ- 
ent) approach was presented by us in Ref. PD| , in which MEC and two-body seagull effects 
were evaluated using two methods: through a procedure of the minimal substitution in the 
NN potential and through a direct diagrammatic computation in the framework of the Bonn 
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OBE picture of the NN interactions. 

It should be said that despite a resemblance of many physical ingredients of Refs. |3"7|-IU 



the results of these works are sometimes rather different, what may indicate unnoticed com- 
putational errors or unjustified approximations. For instance, the differential cross section 
at the forward angle and the photon energy u = 100 MeV found in Ref. [33] is only 2/3 of 



that found in Refs. |38|j4(| (in this comparison, polarizabilities of the nucleon are omitted). 
There are large discrepancies at the backward angle either. Very recently, one more calcu- 
lation in the approach close to that of Ref. |3j| was presented . Their results differ from 



our and previous predictions too, especially at energies to ~ 100 MeV. Possible reasons for 
that are discussed in Section |V| below. 

Recently methods of effective field theories have been applied to deuteron Compton 
scattering as well |42"|-fl4|| . In such an approach, a model- independent part of the low-energy 
scattering amplitude which dominates in the chiral limit of m % — > (but still ^> MA, 
where A = 2.2246 MeV is the deuteron binding energy and M is the nucleon mass) was 
found in a closed analytical form [|42T| . Generally, the results of both calculations are similar 
to those obtained by virtue of the "standard-nuclear-theory" technique. The advantage of 
the calculations |^~H is t na t, they naturally include important non-static effects in the pion 
propagation which is an outside feature for the "standard" theory based on the notion of 
the NN potential. A disadvantage also exists which is related with unavoidable truncation 
of the expansion series leading to a lost of contributions important for a determination of 



the neutron polarizabilities. The A-isobar is one example. See Section [V]] for a further 
discussion. 

In the present paper we complete the calculation with the nonrelativistic Bonn OBE 
NN potential (OBEPR) started earlier [J4U|J . Technically, this is done in the framework of 
the diagrammatic approach which relies on an explicit consideration of relevant Feynman 
diagrams of the reaction in the momentum representation. It avoids Siegert-like transfor- 
mations and it is rather convenient for incorporating non-static and relativistic corrections 



45| , |46[| . Because of inherent restrictions of the potential picture, our analysis covers energies 
below pion photoproduction threshold. 

Our model is essentially nonrelativistic. However, we include a few most important 
relativistic corrections (like the spin-orbit interaction) into the one-body electromagnetic 
current and seagull. After a brief description of the notion of the seagull given in the 
next Section, we introduce the Hamiltonian of the model and analyze one- and two-body 
electromagnetic operators. Then we calculate the Compton scattering amplitude and discuss 
the obtained results. 



II. HAMILTONIANS, CURRENTS AND SEAGULLS 



In the framework of the time-ordered perturbation theory, a computation of the photon 
scattering amplitude starts with a specification of system's effective degrees of freedom and 
the system's Hamiltonian including its dependence on the external electromagnetic 

vector potential A^. We need both linear and quadratic terms in the expansion of H[A] 
in powers of A^ which determine operators of the electromagnetic current j^(x) and the 
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electromagnetic seagull S flv (x, y) for the system and, correspondingly, the so-called resonance 
R and seagull S parts of the Compton scattering amplitude. Simplifying a real situation, 
we write 



H[A](t)=H(t) + J [j,(x)A»(x 
1 

~ 2 



d x 

X0=t 

S, u (x,y)A»(x)A»(y) 



d 3 x d 3 y + . . . 



xo=yo=t 



(2.1) 



where S^ u (x,y) = S Ufl (y,x) is assumed to be a symmetric function of its arguments. Ac- 
cordingly, the photon scattering amplitude of 



\ l ) +1 (e,k)^\f)+ 1 '(e',k') 



(2.2) 



reads 



T(u,8) = R(u,6) + S(u,8) 



to leading order 0(e 2 ), where 



r(uj,6)=y: 



(f\e'^J,(-k')\n)(n\^Uk)\t) 
E„ — Ea — uj — iO 



(e <-> e'*, k «-> -A;') 



(2.3) 



(2.4) 



and 



5(w,0) = (/|^6"S^(-A: , ,A;)|i). 



(2.5) 



Here = k is the photon energy, 9 is the scattering angle, E n are energies of eigen states 
\n) of the system, j^{k) means a Fourier component of the current density, i.e. 



-ik-x 



x =0 



and 



Sfiu{ k , k) 



S^v(x,y) e 



ik' -x—ik-y 



a;o=yo=0 



d 3 x d 3 y. 



(2.6) 



(2.7) 



In a more general situation, the Hamiltonian ([2.1|) can depend on time derivatives of the 
vector potential too (e.g. owing to a presence of terms dependent on the electric field). 
Nothing changes then in Eqs. (|2.4j) - (|2.5|) with the except that the Fourier components of 
the current and seagull become dependent on both the space and time components of the 
photon momenta. 

As is well known, the structure of the effective Hamiltonian (and thus that of the current 
and seagull too) is closely related with a choice of the effective degrees of freedom. In 
the present context of low-energy deuteron Compton scattering, we consider nonrelativistic 
nucleons as the only dynamical variables of the system, whereas all mesons, anti-nucleons 
and other degrees of freedom are encoded into the internal structure and effective interactions 
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of the nucleons themselves. Such an approach is certainly applicable at energies below pion 
threshold. 

With this choice, the resonance amplitude R comes from low-lying (two-nucleon) in- 
termediate excitations n of the deuteron, including the deuteron itself. It corresponds to 
two-step scattering via photon absorption followed by photon emission and vice versa. This 
piece can have and generally has the imaginary part. Meanwhile, the seagull amplitude S is 
real and corresponds to processes, in which the photon absorption and emission happens at 
indistinguishable time moments, as seen at the considered energy scale. Among such pro- 
cesses are excitations of heavier intermediate states like nNN, which describe, in particular, 
meson exchanges between photon interaction points and an internal structure of the nucleon 
related with intermediate-meson production. Considering meson-exchange processes in gen- 
eral as instantaneous (unretarded) and neglecting dependence of meson propagators on the 
photon energy, we retain a retardation correction for the pion exchange which is known 
to be quite important in the seagull amplitude S and to lead to effective modifications of 
nucleon polarizabilities in nuclei |47fl . As for contributions to S due to the nucleon internal 



structure, it is taken into account by introducing polarizabilities of the nucleon. 

Both the current operator j^{x) and the seagull operator S^x, y) have to be consistent 
with the nuclear Hamiltonian H(t) of the system and satisfy conditions of the gauge invari- 
ance. Generally, these conditions take the form of the conservation of the electromagnetic 
current j'^LA^x) found in the presence of the external vector potential A^. 

= PMAKx) = d k Jk [A}(x) +z[H[A}(t), 3o [A}(x)} + d° A j [A](x). (2.8) 

Here the Latin index k runs over the space components and the time derivative d\ acts only 
on the external potential A^. In the simplest case of the time- local Hamiltonian (|2.1|) the 
current j M L4](a;) is given by the three-dimensional variational derivative of -ff[A], 

jM]{x) = 5 ~ImW = Ux) ~ I s ^ y)A ^ y) d3y + • • • ■ (2 - 9) 

In this case the term with d° A in Eq. ( |2.8j ) must vanish, since this is the only piece which de- 
pends on the time derivative of A^. Therefore, the charge density j [^4](^) is A-independent, 
and the following consistency equations emerge 



\jo(x),H(t)) = -i?^ (2.10) 

and 

\jo(x)J*(v))=i dS *^' V \ S M (x,y) = S O0 {x,y) = 0. (2.11) 

Here all operators are taken at the same time moment t = xq = yo- 

In the nonrelativistic approximation, which will be used in the following consideration 
of two-body effects, the charge density jo(x) is not affected by meson exchanges (Siegert's 



7 



theorem) and therefore coincides with the one-body charge density of the two nucleons 
% = 1,2 : 



.70 (*) = Jo 1] ( x ) = H eZi5(x - ri 



l + rf 



(2.12) 



i=l,2 



Then Eqs. ( [2.10| ) and ( |2.11| ) give relations between the nuclear (two-body) potential V 
standing in the nuclear Hamiltonian V and the two-body parts of the current and seagull, 
\ 2] and SE 



J 



and 



'£ ] (x)M =-zV-j[ 2 l(x), 



i 1] (x),jT J (y) 



(2.13) 



(2.14) 



The resonance and seagull contributions at zero energy are constrained by the Thirring's 
low-energy theorem. Within the nonrelativistic accuracy we have 



T(0, 9) = R(0, 6) + S [1] (0, 6) + S [2] (0, 6) 



2„2 



Z 2 e 
AM 



e ■ e 



(2.15) 



where M is the nucleon mass, Z = 1 and A = 2 for the deuteron, and the radiation gauge 



k - e = 0, 



k' • e' = 



(2.16) 



is assumed for the photon polarization vectors. In the absence of the two-body currents 
and seagulls, the model-independent relation (|2.15|) is fulfilled due to a balance between the 
one-body seagull contribution, 



S [1] (0,#) 



M 



e ■ e 



(2.17) 



and the resonance amplitude R(0,6) = (NZe 2 /AM)e ■ e'*. Here N = A — Z. The pres- 
ence of the two-body currents results in an enhancement of the resonance amplitude, viz. 
R(0,6) — ► (NZe 2 /AM)(1 + K)e ■ e'* for spinless nuclei, where k is the enhancement pa- 
rameter standing in the modified Thomas- Reiche-Kuhn sum rule (see, e.g., the review 
for a discussion and further references). Then, in order to support the balance sug- 
gested by the low-energy theorem (|2.15|) , a two-body seagull contribution is required. It is 
S^(0,9) = —(NZe 2 1 AM)ne ■ e'* for the spinless nucleus. For a general case of a nucleus 
of spin S > 1, the two-body seagull amplitude is characterized by the scalar and tensor 
enhancement parameters, k and Kt- 



S [2] {0,6) 



N Z 2 
^ e 2 e + k t SiSj + SjSi - -S(S + 1)% } 



AM 



(2.18) 



Now we proceed with a consideration of free and interacting nucleons. 
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III. HAMILTON! AN FOR A SINGLE POL ARIZ ABLE NUCLEON 



A. Leading-order effects 



Phenomenologically, the dipole polarizabilities a and (3 are defined as low-energy param- 
eters determining the quadratic- in-the- field energy shift Vp Q i, Eq. ( |1.4j ). This shift has to be 
added to a "bare" Hamiltonian i/o|y4] which is linear in the electromagnetic field, describes 
an "unpolarizable" nucleon with the electric charge eZ and anomalous magnetic moment 
en/2M and produces the so-called Born contribution to the Compton scattering amplitude. 
In the relativistic phenomenology, the standard choice for -ffo[^4] and hence the standard 
definition of the unpolarizable nucleon is given by the Dirac-Pauli Hamiltonian 

H Q [A] = eZA Q + oc ■ (p - eZA) + (3M + ^-Bo^ (3.1) 



with F, 



Actually, the given form of i?o[^4] is valid only for the nucleon 
interacting with real photons. This is all we need in the present paper. In the more general 
case, additional derivatives of the electromagnetic field appear in i?o[^4] as wen 0- They 
account for electromagnetic form factors of the nucleon, i.e. its finite size. 

Polarizabilities manifest themselves in low-energy Compton scattering as an 0(uj 2 ) ad- 
dition to the Born amplitude, the latter becoming the Thomson scattering amplitude 
— (e 2 Z 2 /M)e ■ e'* at zero energy. In order to correctly identify the contribution of the 
polarizabilities, O(oj 2 ) terms in the Born amplitude have to be retained as well. Since some 
of them are of order 0(uj 2 /M 3 ), an effective low-energy Hamiltonian covering all the 0(uj 2 ) 
terms has to include relativistic corrections up to order 0(M~ Z ). 

A nonrelativistic reduction of the Dirac-Pauli Hamiltonian valid to order required was 
found in Ref. 



IE 



Using the Foldy-Wouthuysen method |50| , |5lf or expelling lower com- 

a lengthy but straightforward 



ponents and higher derivatives as described in Refs 
computation gives^ 



H [A] = eZA + 



7V- 



7T 



2\2 



e(Z + n) 



2M 8M 3 
e(Z + 2n 



a- H 



8M 2 



2M 

V-E + cr-(Ex7r-7rxE) 



+ -^-{7r 2 , cr ■ H} + -^r{cr ■ ir.n ■ H} 

8M 3 1 ' 1 8M 3 1 ' 1 



+ 



+ 



16M 3 

„2 



8M 3 



{tt, V x H - E} + (<r x V) • (V x H - E) 

Z 2 + Zk + k 2 )~E 2 - Z 2 B. 2 } + 0(M- 4 ). 



(3.2) 



We give the answer in the form obtained in Ref. IS]. In Ref. [p(j] the anomalous magnetic 
moment is not considered and the final result contains a sign mistake. In Ref. terms of order 
0(e n ), n < 2 were only retained. 
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Here 7r = p — eZA is a covariant momentum, {A, B} denotes the symmetrized product 
AB + BA, and E means the time derivative of the electric field. In the region lying outside 
any sources J p of the electromagnetic field, the combinations V • E = J and V x H — E = J 
vanish, so that the above equation turns out even simpler. 

When anti-nucleon degrees of freedom are removed and absorbed into new effective inter- 
actions, the resulting effective Hamiltonian ( |3.2|) becomes non-linear in the electromagnetic 
field. In particular, it contains polarizability-like parts which have to be kept in compu- 
tations using nonrelativistic variables alone. Among these parts is the term k 2 E 2 which 
imitates a negative electric polarizability of the neutron and which is known due to Foldy 



One can easily check that the Hamiltonian ( |3.2| ) exactly reproduces the Born amplitude of 
nucleon Compton scattering to order 0(u 2 ) which is explicitly given, e.g., in Ref. ||. More- 
over, all the 0(uj 2 ) terms in the scattering amplitude are retained when the kinetic energy 
in the nucleon propagator is calculated to leading order 0(M~ l ) (i.e. nonrelativistically), 
the electromagnetic current is taken to order 0(M~ 2 ) (i.e. with the spin-orbit correction), 
and the full electromagnetic seagull up to order 0(M~ 3 ) is taken as it stands in Eq. (|3.2|) . 

In the present work we adopt a few further simplifications to Eq. ( |3.2j ). First, we ne- 
glect those parts of the Hamiltonian which do not contribute to the 0(u 2 ) terms at all. 
These are the 0(M~ 3 ) parts of the kinetic energy and the current. Second, we neglect 
the C(M" 3 ) part of the spin-dependent seagull which gives an 0{u 2 ) contribution to the 
Compton scattering amplitude but does not contribute to the differential cross section of 
Compton scattering to order 0{uj 2 ) with unpolarized nucleons. Third, omitting the 0(M~ 3 ) 
component of the kinetic energy, we omit also a ~ e 2 Z 2 /M 3 part of the seagull standing; in 
— 7r 4 /8M 3 ; moreover, we omit such a part in the coefficients of the fields squared. Thus, 
we use the following effective Hamiltonian for a single nucleon which interacts with real 
photons: 

HW[A] = eZA + 4 ~ <Z + , K) * ■ H ~ <Z + , 2/6) <r • (Ex.-.xE) 
L J 2M 2M 8M 2 v ; 

- -An(a + 5a )E 2 - -4tt/3H 2 , (3.3) 
2 2 

wheref] 

e 2 k 2 + Zk _ J —0.85, proton, , . 

da ° = 4M 3 ~ I -0.62, neutron. ^ ' 



4 The correction <5«o has a direct relation with the difference mentioned in Section | between a n 
and the "static polarizability" a n found in the neutron transmission experiments [2C ,21] . In fact, 
in the formalism used in these works a n denotes the coefficient of E 2 in the effective nonrelativistic 
Hamiltonian. In order to get warning against wrong generalizations note, however, that the coeffi- 
cient of E 2 in the proton case is not equal to the static polarizability a p which differs from a p by 
a term containing the electric radius of the proton 
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The corresponding electromagnetic vertices, i.e. matrix elements of the one-body current 
and seagull in the momentum representation, read 

efj W (k; p',p) = -|^e • (p + p') - (Z + k) iua • s 

[Z + 2k) iuja ■ e x (p + p') (3.5) 



and 



1M 2 



s'*^S$(-k', k) = e .e'* + ^(Z+ 2k){u + J) ia 



+ 47iW(a + 5a ) e ■ e* + Attlulu'/3 s ■ s'*, (3.6) 

where u and u' are the initial and final photon energies, 

s = k x e, s' = k x e' (3.7) 

are the magnetic field vectors, and we have used the radiation gauge ( j2.16 ) for the photon 
polarization vectors. It is worth mentioning that the absence of the kinetic term —p 4 /8M 3 
in the Hamiltonian ( |3.3| ) allows us to use self-consistently nonrelativistic phenomenological 
potentials developed for a description of iViV interactions at low energies. 

The Hamiltonian ( |3.3| ) with the leading relativistic corrections included possesses an 
accuracy of about e 2 /167rM 3 = 0.17 x 10~ 4 fm 3 for treating the leading-order effects of the 
polarizabilities. For example, being used in the lab frame, the Hamiltonian ifM[y4] and the 
vertices (|3.5|) - (|3.6|) generate the following scattering amplitude at the forward angle, 

T lab (£ = 0°) = ( - ^ + 4™ 2 (a + P))e- e'* 

2 2 

r ~ K xua-e!* x e + C(cj 3 ). (3.8) 



2M 2 

The 5ao and spin-orbit contributions of the seagull properly correct w-dependent terms 
coming from the resonance amplitude R and bring the resulting amplitude (|3.8| ) into an 
exact agreement with a known low-energy expansion of T (given, e.g., in Ref. ||25|). At the 
backward angle, the amplitude found with the Hamiltonian (|3.3|) reads 

T Uh {6 = 180°) = ( - ^ + 4W(a - $)) e • e'* 

e 2 U 2 + 4Zk + 2Z 2 ) , „ M ^, . . 

+ AM 2 Z{UJ + J)CT ' 6 6 + ° {UJ ) - (3 ' 9) 

This time an exact result is slightly different. It contains an additional term 
(e 2 Z 2 /2M 3 ) ujuj'e-e 1 * which comes from a recoil (9(u; 2 )-correction to the Thomson amplitude 
and which is lost in Eq. ( |3.9|) because of omitting the e 2 Z 2 /M 3 pieces of the seagull. | 



5 Making such a comparison of the two amplitudes, one has to take into account a different 
normalization of the nucleon states. It is one particle per unit volume in the present paper and 2E 
particles in Ref. |2j| • 
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B. Polarizabilities and the Baldin sum rule 



In the case of scattering, the seagull amplitudes Q2.5|) for the proton and neutron 
contribute coherently and dominate the scattering amplitude (|2.3|) at energies of a few tens 
MeV. Their joint result depends only on the isospin-averaged polarizabilities of the nucleon, 
viz. a N = \(a p + a n ) and f3 N = \{f3 p + f3 n ). 

In the following we will consider the difference — (3n as the only free parameter of the 
nucleon structure. It is hard to reliably predict this difference, because it can be affected by 
t-channel exchanges with poorly known couplings (like the u-meson exchange) - see, e.g., 
Refs. [0,0. Meanwhile the sum a^ + ^N can De safely found from the well- convergent Baldin 
sum rule ( |1.3|) . This is quite sufficient in the present context. 

There are several evaluations of the dispersion integral in Eq. ( |1.3| ). Earlier calculations 
gave a p + f3 p = 14.2± (0.3-0.5) lj5|-|6l and a n + (3 n = 15.8 ± 0.5 @ (we comment on the 
other result, a n + j3 n = 13.3 []5B| below). A recent re-analysis |S7[] gave lower values: 



a p + (3 p = 13.69 ±0.14, (3.10a) 
a n + (3 n = 14.40 ± 0.66. (3. 10b) 

Doing our own calculations with modern sets of photoabsorption data, we also obtain some- 
what lower values than those found in 70's. However, they are not so low as those in Ref. 
[57j, especially for the neutron. 

Specifically, finding cr tot (u;) through the set of pion photoproduction amplitudes of Ref. 



5q| at energies below 400 MeV, taking total photoabsorption cross sections from Refs. |59lj6 



at energies 0.5-1.5 GeV, making a smooth mixture of the "theoretical" [^8] and experimental 



[ 59| , |60|| cross sections in between, and using a Regge parameterization of a tot (u)) at energies 



oj > 1.5 GeV (the same as in Refs. [|T],|55[), we obtain 

a p + (3 p = U.0, (3.11a) 
a n + p n = 15.2. (3.11b) 

Uncertainties in these numbers mainly originate from the region of u ^ 500 MeV which 
essentially saturates the dispersion integral. They can be again conservatively estimated as 
±(0.3—0.5). For example, we obtain very close results (13.8 and 15.2, respectively) using in 
this computation photo-pion amplitudes from the code SAID |61[ (as of beginning of 1999) 
instead of the amplitudes from Ref. |58|| . The lower value for a p ± (3 pi which follows from the 



SAID amplitudes, is mainly caused by that the pion photoproduction multipole E 0+ (ir + n), 
as given by that partial-wave analysis close to pion threshold, is by ~ 12% too low |62| in 
comparison with predictions of independent analyses like p3[ and with predictions of chiral 
perturbation theory |63|.Q In accordance with (|3.11| ), we accept the following number for 



the isospin-averaged sum of the dipole polarizabilities of the nucleon: 



6 Actually, the SAID authors gave an explicit warning against using the SAID amplitudes very 



close to pion threshold [61|. Addition written after the first submission of the present paper: Very 
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ol 7v + Pn = 14.6 



(3.12) 



There are several reasons why we prefer to rely on our own findings ( |3.11| ) both for the 
proton and the neutron rather than on the quoted recent results fl3.10| ). For the proton case, 
the central number for the sum of the polarizabilities obtained by the authors of Ref. p? 



is 



shifted down by their use of the SAID amplitudes very close to pion threshold. This shift 
almost explains the difference between ( |3.10a| ) and ( |3.11a| ). It is worth saying that the tiny 



uncertainty ±0.14 ascribed to a p + (3 P in Eq. Q3.10a| ) represents only statistical errors in the 
experimental data on cr tot . It does not include systematic errors which are equal to 2—3% 
in crtot and hence produce the uncertainty of at least ±0.3 in a p + /3 P . 

For the neutron, we are even more far from reproducing the very low value obtained in 
Ref. 1 15 711; we are also far from the result of Ref. 15611, where the number obtained was even 



lower. The reason might be in a different use of the (indirect) data |59[ on the neutron total 
photoabsorption cross section cr to t(7^)- Close to the A-resonance energy, the cross section 
cr tot(7 ri ) given in Ref. |59| is by ~ 20% (!) lower than predictions of all modern partial- wave 
analyses of pion photoproduction. The procedure used in Ref. [[59] to extract o" to t(7 n ) from 
the primary cross section <r tot (7(i) obtained with the deuteron target is not so clear in the 
A-resonance region, in which medium corrections are large. That is why we believe that the 
results of Ref. |59 for the neutron should not be taken seriously at energies below 400 MeV. 
As was already said, in our own evaluation of Eq. ( |3.11b| ) we have found o" to t(7^) below 400 
MeV through the partial-wave analyses [|58|,pl[ . 



C. Higher-order corrections 



It is clear that higher-order kinematical corrections neglected in Eq. ( |3.3j ) are suppressed 
by powers of u/M and therefore they are small below pion threshold.^] An actual accuracy 
of the effective Hamiltonian ( |3.3| ) is determined by dynamical effects which originate from 
the pion and A-isobar structure of the nucleon and give corrections of the relative order 
(uj/m-n) 2 . They become important at energies ^ 70 MeV. To next-to-leading order, such 
effects are parametrized by eight structure constants of the nucleon, viz. the quadrupole 
(atE2, Pm2), dispersion (a Eu , fi\iv) and spin (j E1 , j m , ^y E2 , j M 2) polarizabilities of the 
nucleon |64]- p6|j25| , as represented by the following effective Hamiltonian [p5|] : 

5H®[A] = -Un(a E „E 2 + (3 M ^) - ^{a m E% + PmiB.%) 
— -47r(7£icr ■ E x E ± Jmi& HxH 



recently a new set of the SAID amplitudes appeared on the SAID web-site (solution SM99K). This 
set is free from most above-mentioned near-threshold problems, and its use in the evaluation of 
the Baldin sum rule gives results which perfectly agree with the estimates (3.11). 

7 The suppression in T-even observables like the differential cross section is in fact (to/M) 2 . 
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2'j E2 E ij a i H j + 2^ M2 H ij a i EA. 



(3.13) 



Here 



E ij = oi^iEj + VjEi), Hij = -(ViHj + VjHi 



(3.14) 



are quadrupole strengths of the electric and magnetic fields. Such an effective interaction 
contributes to the seagull amplitude of 7iV scattering which gets an addition 

e'*^5S [ ^(-k', k) = 47iW[e • e'*5a + s ■ s'*5(3 
i 



(3.15) 



+ - (u + J) (a- ■ e* x e (~/ M2 - j E i) + cr ■ s'* x s (j E2 - 7mi 
— i(cr ■ ke • s'* — cr ■ k' s ■ e'*) r )E2 — *(c • k' e ■ s'* — cr ■ ks ■ e'*)~/M2 



The functions 5a, 5(3 which depend on the photon energies and on the cosine z = k ■ k of 
the scattering angle, 



5a = uu'(a Eu + ^a E2 - 

5(3 = Uu'(j3 M v + ~ TK a E2) 



(3.16) 



can be handled as dynamical corrections to the dipole polarizabilities a, (3 standing in Eq. 

(PD- 

Using estimates for the isospin-averaged polarizabilities of the nucleon obtained through 
fixed-t dispersion relations [25] , 

{a Ev ) N ~ -3.1, {p Mv ) N ^ 9.1, (^2^^27.3, {(3 M2 )n ^ -23.0 (3.17) 

(units are 10 -4 fm 5 ), we find, e.g., that the contribution of 5a, 5(3 increases the backward- 
angle differential cross section of 7<i scattering and makes the same change as a shift of 
a N - (3n by -1, -2 and -4 at 50 MeV, 70 MeV and 100 MeV, respectively. 

It is not easy to estimate a model uncertainty in the above numbers Q3.17| ). In particular, 
they are sensitive to the so-called asymptotic contribution Af s to the Compton scattering 
amplitude A\ which was represented by the a-exchange of the effective mass m a = 500—600 
MeV in the framework of Refs. |25],29|. Recently, an alternative dispersion approach was 
presented |30| , |6"7| which allowed to avoid an explicit introduction of the a-exchange and 
to calculate the amplitude A\ and the quadrupole polarizabilities of the nucleon under 
reasonable assumptions on the reactions 77 — > tttt and tttt — > NN!^ Specific numbers 
obtained in Ref. |)7] for the higher-order polarizabilities (given for the proton case only) 



8 There are many predecessors of this approach, of which we would like to indicate Refs. [2^.27,68] 
as most recent works, in which further references can be found. 
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are rather close to those obtained earlier f25f . Their use for evaluating the deviation of 
the backward Compton scattering amplitude from the low-energy expansion of order 0(u 2 ), 
Eq. (|3.9|), gives only a 3% bigger effect than that obtained with our numbers ( |3.17| ). More 
cautiously, we could state that the effect of the quadrupole and dispersion polarizabilities 
is known within 20%, where the last number is obtained by a reasonable variation of the 
effective a-meson mass. 

In order to evaluate the spin-dependent contribution in Eq. ( |3.15| ), we use spin polariz- 
abilities of the nucleon as found through the dispersion relations too [^,^,^] : 



(iei)n ^ -3.7, (7mi)jv ^ 2.3, {^ E 2)n ^ 1-4, (jm2)n ^ 0.6 (3.18) 



(units are 10 4 fm 4 ). Writing Eq. ( |3.18 ), we have corrected predictions for 7's of Refs. 



which include a poorly constrained asymptotic contribution Af arising in the fixed-t 
dispersion relation for the invariant amplitude A 2 of nucleon Compton scattering. Since 
A™ determines the backward spin polarizability of the nucleon, j n = —^ei + 7mi + lE2 — 
7M2, which was recently re-evaluated through a more reliable backward dispersion relation 
we have introduced the appropriate changes to 7's. Specifically, they are 5(jei)n — 
-8{jmi)n = -8{ r YE2)N = S{1m2)n = -\Ki*)n, where 5{^) N ~ -4 is a correction to the 
previous estimate |]25|j62]] of (7^)^- About one half of that corrections stems from the 77 and 
7/ exchanges. At backward angles, the spin effects of order 0(u 3 ) result in enhancing the 
coefficient (e 2 /4M 2 )(/t 2 + AZk + 2Z 2 ) in Eq. (U) by 2iruu% with (7^ ^ 7 and make an 
increase in the differential cross section of backward-angle 7c? scattering which is about one 
third of what the (5a, 5(3) correction does. 

Recently, there was a controversy on the value of j n , at least for the proton. It was 
experimentally found |13[ that (j^p = — 27.1±2.2 jj^ (i n units of 10~ 4 fm 4 ) and, therefore, it 



largely deviates from theoretical predictions [F2^,^,^,[T{J which give ('j^p = —36.7 to —39.5 
and (7, r )„ = 50.3 to 52.5. Taking this deviation seriously and assuming that the isospin- 
averaged backward spin polarizability (7^)^ may differ from the above-accepted value of 
(Jtt)n = 7 by as much as +10, we should conclude that the theoretical differential cross 
section of •yd scattering at backward angles is visibly enhanced due to the effect of (j^n- 
Accordingly, the value of 6ln— Pn extracted, for instance, from 100 MeV data on 70? scattering 
should be shifted up by as much as +2 x 10 -4 fm 3 , what is not negligible! On the other hand, 
recent Mainz measurements of backward-angle proton Compton scattering J7T| done with 



the deuterium target give the differential cross section which is lower than that obtained at 
LEGS JI3"| , and these newer cross sections assume that the theoretical value of (7^ ~ —37 
is fully compatible with the data. Moreover, recent Mainz measurements of proton Compton 



scattering data done with the hydrogen target |F2 seem to fully exclude the previous finding 



fi"3| as well and to completely agree with the quoted theoretical predictions. In view of that, 
we rely our following analysis on the theoretical values of the spin polarizabilities of the 
nucleon specified by Eq. ( |3.18| ) and assume that the related theoretical uncertainties in 
extracting — Pn are less than w 0.5 x 10 -4 fm 3 . 

At forward angles, all the higher-order corrections ( 3.1 5|) are less important. 
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IV. TWO-BODY CURRENTS AND SEAGULLS 



A. Potential-induced electromagnetic currents and seagulls 



The remaining part of the Hamiltonian H[A] is related with two-body interactions of the 
nucleons. In the absence of the electromagnetic fields, such interactions can be represented 
by a (generally non-local) iViV-potential V which has to accurately describe differential cross 
sections and polarization observables in NN scattering at energies below pion threshold, 
as well as the deuteron binding energy. There are several phenomenological potentials of 
that sort in the literature. We have chosen to use the Bonn potential (specifically, its 
nonrelativistic version OBEPR) |73,73J], because it implies a very simple physical picture of 



one-boson exchanges (OBE) mediating the NN interaction and, in the framework of this 
picture, allows constructing the meson-exchange current jffl and the meson-exchange seagull 
SjfJ directly from the corresponding Feynman diagrams. Of course, the OBE picture cannot 
be true in all detail. However, at least, it takes fully into account the most important 
long-range contribution, which is the one-pion exchange. 

In the momentum representation, the OBEPR potential has the form 

v(p[, p' 2 ; Pl , Pa ) = £ v«( P 'i, p' 2 ; Pi, Pa), (4-1) 

a=7r,?j,(5,CT,cj,p 

where and p' { are the initial and final momenta of the i-th nucleon subject to the constraint 
Pi + P2 = Pi + P25 an d V a are potentials stemming from the exchanges with the specified 
mesons a — 7T, rj, 6 (which is a (980) in the modern notation), a, u and p (see Fig. fflo). 




FIG. 1. The one-boson-exchange potential (a), and the corresponding meson-exchange current 
(b) and meson-exchange seagull (c). Strong and radiative aNN vertices include antinucleons (due 
to PS couplings and/or relativistic corrections) and form factors. 

Let us consider in some details the pion exchange which determines the long-range part 
of V and gives the biggest contribution to the matrix elements of MEC and MES relevant 
to low-energy 7c? scattering. The potential V 11 is velocity- independent, i.e. it depends only 
on the momentum transfer q: 

V n (cO = -JL CTl . qo - 2 . qTl . T2 GV(q). (4.2) 
Here g w is the 7rNN coupling constant, and the function G n , 
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GV(q) 



^ 2 (q) _ 

q 2 + m 2 



(4.3) 



contains the pion propagator and the ttNN vertex form factor of the monopole form, 

pn A 2 — ml 



A 2 +q 2 



(4.4) 



given by the isospin- averaged mass of the pion, m n , and the cutoff parameter A n . 

The pion-exchange current j* can be obtained by attaching the photon line to the ex- 
change pion and to the 7rNN vertices as shown in Fig. in which the electromagnetic 
^/7tNN vertex arises due to a momentum dependence of the ttNN coupling. This mo- 
mentum dependence comes partly from the derivative, or the factor of q, standing in the 
nNN vertex (or, equivalently, from a contribution of antinucleons in the formalism of the 
pseudo-scalar coupling adopted in Refs. |73| , |74]| ). Then the minimal substitution 



V(r ■ 7r) — ► V(r ■ 7r) - ieA 



T ■ 7T 



(4.5) 



in the effective ttNN Lagrangian generates the well-known Kroll-Ruderman component of 
the jirNN vertex. An additional momentum dependence is introduced by the vertex form 
factor, Eq. ( (4.4|) , and it should also be taken into account. 

Without knowing the dynamical nature of F n (q), there is no unique way to restore the 
electromagnetic current associated with the form factor. Different prescriptions proposed in 
the literature for maintaining gauge invariance in such cases (see, e.g., Refs. |75|-|83|) give 
different answers, especially in the region of high momenta q. Fortunately, at low momenta 
g< A x which are only relevant to the present consideration, such ambiguities are expected 
to be small, as is suggested by the Siegert's theorem. In the following we choose a simple 
phenomenological way explicitly formulated by Riska [77| and Mathiot |78||.fl That is, we 
assume that the vertex form factor ( f4.4|) results from a propagation of a fictitious particle 
A of the mass A„-, which has the same quantum numbers as the pion and mediates the 
pion interaction with the nucleon (see Figs. |2]a and |p). Accordingly, the photon interacts 
with the particle A as well (Fig. |2|&), and this gives an additional contribution to the vertex 



'ynNN which restores the fulfillment of the generalized Ward-Takahashi identities [[84]] for 
the transition amplitude of 7 A — > ttN and restores the electromagnetic current conservation 
in the meson-nucleon system. 



9 In essence, the solution given by Riska and Mathiot was earlier obtained by Arenhovel [7q| , who 
used the minimal substitution (4.5) and considered the specific case when the monopole form factor 
( |4-4| ) appears to the first power in the potential V" K (q). A straightforward generalization to the case 
when this form factor is squared, as in Eq. (|4.3| ), is easily derived through a differentiation with 
respect to |83[| . After this differentiation, the Arenhovel's prescription [75 becomes identical 
to that proposed by Riska and Mathiot Wl\ffl 
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a 



i i 
l I 
I I 



a) 



b) 



Xx 



+ crossed 



FIG. 2. Diagram a: the meson- nucleoli form factor F a (c[) viewed as resulting from an intermedi- 
ate heavy-boson exchange A. Diagrams b and c: corresponding contributions to the electromagnetic 
meson- nucleon vertices. 



Evaluating the diagrams shown in Fig. |l|fe with the vertices shown in Fig. 06 and with 
the static propagators of all particles,P] we obtain the well-known result for the pion MEC: 



j^(k;p'i,P2;Pi,P 2 ) = -ie(Ti x r 2 ) 



AM 2 

+ ie(qi - q 2 ) (ti x t 2 



o-i (0-2 ■ q 2 ) G r 7 r(q 2 ) -(1^2) 



4M 2 



<T! • • q 2 Gi u (q 1 ,q 2 ) 



(4.6) 



Here we introduced the function (cf. Refs. |]77|j78| , |83| 

^(qi)^(q 2 ) 



Ci7r(qi,q 2 ) 



(q? + ml) (q| + ml] 



+ ql + Al + q? + A2 



(4.7) 



which provides a combination of propagators of the exchanged pion and the particle A as 
they appear in the case of a line with one electromagnetic vertex (see Fig. |3]&). The vectors 



qi = pi - Pi, q 2 = p' 2 - p 2 



(4.8) 



are the momenta transferred to the nucleons. These momenta are subject to the constraint 
q x + q 2 = k, where k is the momentum of the incoming photon. Using the identity 



GU<li) - GV(q 2 ) = (q 2 - q?) G\^i, q 2 ) 



(4.9) 



one can easily check that the pion-exchange current ft4.6|) satisfies Eq. Q2.13|) , which has the 
following form in the momentum representation: 



k • ,T(k; p[, p 2 ; p 1} p 2 ) = eV n (p[, p 2 ; p 1 + k, p 2 )Z 1 

- e Z 1 ^(p , 1 -k,p 2 ; Pl ,p 2 ) + (l 



2). 



(4.10) 



Here Zi (i = 1, 2) are the electric charges of the first and second nucleon, Eq. fl2.12|) . 



10 Non-static, i.e. retardation, corrections will be considered in the next subsection. 
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G a (q) G la (^ u q 2 ) 

4 1 I ]* 



a) 



d i h 



G 2a (qi,q 2 ,K) 

in i j f 

I I I VI \i I 



FIG. 3. Modifications of a light-boson propagator by intermediate heavy bosons A. Shown are 
the cases with zero, one, and two photon vertices on the line. 

The pion-exchange seagull S* is determined by the diagrams shown in Fig. |l]c. There, 
the electromagnetic meson-nucleon vertices include again the form factors generated by the 
fictitious particle A through the mechanism shown in Figs. |2]& and ^]c. In the case of the 
pion exchange the very first diagram of Fig. |2|c is absent. However, it contributes when 
the strong meson-nucleon vertex of the OBE potential has a quadratic dependence on the 
particle momenta. This is the case for a, 5, u and p exchanges (see Ref. [74J, Appendix 
A. 3). Evaluating the diagrams |T]c, we obtain 



e 2 9l 



e i* e j^(- k ' ; k ; p'i, P2; Pi> P2; = ± 12 



Tv 



AM 2 



e 2 9l 



T 12 { 



(T l e(T 2 e'* G^Kx) + (1 2) 

(T 1 ■ e cr 2 • q 2 (q 2 - K x ) • e'* G ln (K u q 2 ) + (1 2) 
+ (e^e'*, K!^-K 2 )} 



+ AAP Tl2<Tl ' qi °" 2 ' q2 ^(<ll'<l2> K l> K 2)- 



(4.11) 



Here q x and q 2 are again given by Eq. 

Ki = q x — k 
K 2 = q 2 - k 

The isotopic factor T 12 is equal to 



, and the vectors Ki and K 2 are defined as 



q 2 
qi 



k'. 
k'. 



T 12 — T 21 — Ti ■ r 2 — r^r ; 



z^z 
2 ■ 



(4.12) 



(4.13) 



The function D n is proportional to the amplitude of pion Compton scattering modified by 
the form factor corrections. It reads 



Ar(qi, q 2 , Ki, K 2 ) = 2e • e* G ln (q 1 , q 2 ) 



+ 



( qi + KO • e (q 2 - Ki) ■ e'* GW(q 1; q 2 , Id) + (1 ~ 2) 



(4.14) 



Here the function 
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GW(qi,q 2 ,K) = 



^(qi)^(q 2 ) 

I) (q| + ml) (K 2 + 




q 2 + ml q 2 . + ml \ K 2 + m 2 
^TA^ + "q|+A^ J K 2 + Al 



(4.15) 



provides a combination of propagators of the exchanged pion and the particle A as they 
appear in the case of a line with two electromagnetic vertices (see Fig. [5]c). Writing Eq. 
( PH ), we did not assume any special gauge for the photon polarizations. Therefore, that 
equation specifies all individual components of the tensor Sfj. Using Eq. (|4.9|) and the 
identity 



Glutei, q 2 ) - G lw (K, q 2 ) = (K 2 - q 2 ) G 27r (q 1 , q 2 , K) 



(4.16) 



one can check that the obtained MES satisfies Eq. fl2.14| ). In the momentum representation, 

-k', k; p' x , p' 2 ; p 1; p 2 ) = ejf(-k'; p' 1; p' 2 ; p x + k, p 2 )Z 1 

- eZ 1 jf(-k'; p[ - k, p' 2 ; Pl , p 2 ) + (1 <- 2). (4.17) 



We may note that formulas for the seagull Sfj (in the r-space) were derived long ago in 
Refs. P5] , |T7| by considering the appropriate Feynman diagrams, and in Ref. \W\ by using 
the minimal substitution. Neither of those considerations, however, takes into account the 
ttNN vertex form factor. Therefore, in order to achieve a consistency with the pion-exchange 
potential (fO|), we do need Eq. (|47TT| ). 

and seagulls related with other bosons a of the OBE 



Meson-exchange currents j c 



potential can be derived in a similar way. Formulas for j a were already obtained in Ref. [46 
Newer results for Sfj are given in Appendix |A|. 

All the considered MECs and MESs can be called potential-induced, because they contain 
only those pieces which are intimately related with the OBE potential itself and which are 
needed to fulfil the electromagnetic current conservation in the AW system as given by Eqs. 



1|), (|2.13|) and (|2.14|) . Non-potential contributions to and S\f also exist, and now we 



proceed with a consideration of them. 



B. Non-potential contributions 



The most important degree of freedom explicitly missing in the OBE-potential picture 
of the A^V interaction at low energies is an excitation of the intermediate A-isobar. Never- 
theless, within the purely hadron sector (viz. AW) effects of the A-excitation are indirectly 
included owing to the use of fitted parameters adjusted to the experimental data on AW 
scattering. Then, in accordance with the Siegert's theorem, the electric contributions to 
MEC and MES found with such parameters take the effects of the A into account as well, 
at least in the long wave-length approximation. 

This mechanism, however, does not work for magnetic contributions to MEC and MES 
which have to be added independently. The dominating (long-range) parts of such contribu- 
tions come from the one-pion exchange, and they appear through modifications T^, , 
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of the effective 'yirNN and 'yytrNN vertices caused by the 7./VA transition (see Fig. ^). We 
do not include similar modifications for the case of the p meson, because they are com- 
pletely negligible in the present context. It is worth noticing that the vertex T^ n appears 
due to the momentum dependence of the 7riVA coupling, and it is needed to maintain the 
gauge-invariance of the resulting Compton scattering amplitude. 

■pA -pA 

fa 77a. 

. r. 



- 77 a 



+ crossed 

r 1 11 1 " 1 11 

a) b) 

FIG. 4. The effective jaNN and jjaNN vertices. For a = ir they take into account the 
A-resonance contribution. 

Being used for evaluation of the diagrams in Fig. |1]&, the T^ n component of the jtiNN 
vertex gives the following contribution to MEC: 

j (k; Pl , p 2 ; Pl , p 2 ) = 36 ^ 2 — (2r 2 - z(n x r 2 ) 



M A - M -u 



x (2q 2 — i(Ti x q 2 ) x k + (ri — > — Ti, cr 1 — > —er 1 , uj — > —u) + (1 <-> 2). (4.18) 
Writing this equation, we assumed that the form factor of the ttNA vertex was equal to 



that of the nNN vertex. The mass and couplings of the A are taken to be 86 



M A = 1225 MeV, <£ = Q.282e Ma + M , = —. (4.19) 



Actually, the crossed term in ( 4.181) , i.e. the term having Ma — M + u in the denominator, 



vanishes when the operator j*^ acts upon the deuteron state which has the isospin 1 = 0. 

The A-isobar contributes to the seagull operator too. This happens in two ways. First, 
the component of the jtiNN vertex works in the diagrams with one or two contact 
single-photon vertices shown in Fig. [l]c. This gives a contribution which, together with 
pieces without A, can be written through (off-shell) pion photoproduction amplitudes as 

£ :w-(-^k;p;,p- 2 ;p 1 ,p 2 ) = ™ - - iK) + a » 2). (4.20) 

where Ki is given by Eq. (gig ), and the sum over the pion's isospin index a is assumed (see 
Figs. [S]a and |B|). Second, the Tz~ component of the •yyirNN vertex works in the first two 
diagrams of Fig. [l|c with the contact two-photon vertex. This gives the contribution shown 
in Fig. Ub: 



21 



,* C 7rA( c )/ 1/1 / / \ -rp f<r 2 • q 2 C 7r (q 2 ) 

e t eA, (-k , k; Pl , p 2 ; Pl , p 2 ) = { — _ - - - 

x 2uoe!l ■ s + 2o/e A ■ s'* - icocr l ■ e'l x s + ioj'&i • e A x s'* 

+ (1 «-> 2)} + (e «-> e'*, k «-> -k', to -> -a;'). (4.21) 

Here s and s' are given by Eq. ( |3.7|) . The meson-baryon form factor and the appropriate 
electromagnetic coupling with the fictitious A-boson (Fig. 06) are included into Eq. ( [4.21 ) 
via the use of the vertex function G 7r (q 2 ) and the modified polarization vectors 



e Kx 



e A = e - 2Ki T , 2 , A \ (4.22) 



(and the same for e A ; under the crossing, — > — K 2 and s <-> s'*). 



¥y Y 

«JS = la + _j_ + (1 « 2) + 

a) 6) c) 

FIG. 5. The diagrammatic representation of the two-body seagull S^. The content of the pion 
photoproduction amplitude T 77r is explained in Fig. |6|. 
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-O- = V + — r- 

FIG. 6. The amplitude of pion photoproduction. 

Another feature missing in the potential picture of the AW interaction is that the 
exchange-boson fields are generally non-static. Non-static, or retardation, effects are most 
important for the pion exchange due to its large range. It was recently emphasized by Htitt 
and Milstein in their studies of Compton scattering by heavy nuclei [f|7] that the retardation 
correction gives a noticeable contribution to 7T-MES and to the Compton scattering ampli- 
tude. In the framework of our formalism, we take this correction into account by using the 
retarded pion propagator in Eq. (|4.20|) : 



1 1 1 u 2 . 

(4.23) 



K 2 + m 2 K 2 + m 2 -uj 2 K + ml (K + m 2 ) 2 
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Here we neglect the energy carried by the nucleons and replace the pion energy by the photon 
one, us. The adopted procedure is not fully self-consistent, because we neglect retardation 
corrections to the 7T-MEC and to the pion-exchange potential V n . However, the omitted 
corrections are expected to be less significant than the retardation correction to the seagull 
amplitude (cf. Refs. J47|,|87|). 



V. COMPUTATION OF THE AMPLITUDES AND CROSS SECTIONS 

We do actual computations of the scattering amplitude T(E y , 7 ) in the center-of-mass 
frame of the 7c? system. Accordingly, uo = to' and k, k' mean the energy and momenta of 
the photons in the CM frame; also, 7 means the CM scattering angle. The symbol E 1 is 
reserved to denote the photon beam energy in the Lab frame. To specify polarizations of the 
particles, we introduce the helicities of the photons (viz. A and A') and the spin projections 
of the deuteron and nucleons to the beam direction e 2 (viz. m, m x and m 2 ). In the radiation 
gauge ( gTp , 



1 M ■ N / 1 



(Xe x + ie y ), e'= — -=(\'e x r + ie y ), (5.1) 



where the axis x' is orthogonal to k' and lies in the reaction plane xz. 

Using the seagull operator SiL(— k\ k) specified in Section |TJ, we obtain the one-body 
seagull amplitude S^(E y , © 7 ) through a loop integral in the momentum space (see Fig. |7|a): 

(A>'|SW(£ 7 ,0 7 )|A,m> = / * ^^(p-Itf) *- m2 (p-Ik 

x (m , 1 ,m2\e'*e j S [ i} j {-k',k)\rn 1 ,m 2 ) + (1 <-> 2). (5.2) 

The notation here is that the state-vectors like \mi, m 2 ) are used to designate only spin 
variables of the particles. The momentum variables, if any, are indicated as arguments of 
the operators. The sum is always taken over spin projections of intermediate nucleons. The 
subscript 1 in Sj 1 ^— k', k) points out that this operator acts on the nucleon 1. 

The deuteron wave function m2 (p) depends on the relative momentum p of the 
nucleons. For the nonrelativistic Bonn potential, one can use analytical parameterizations 



of x J / mim 2 (p) §i ven i n R- e f- (HI- Note that the authors of the Bonn potential published 
three nonrelativistic versions of that potential which we label OBEPR(A), OBEPR(B) and 
simply OBEPR. The only difference between these versions is in boson's masses, couplings 
and form factors used. The parameters of OBEPR are specified in Ref. 7q| , Table 14. 



Parameters of OBEPR(A) and OBEPR(B) are given in Ref. 0, Table A.3, part A and 
part B, respectively. We always give our predictions obtained with the OBEPR version 
unless other stated explicitly. 

Similarly to Eq. (|5.2p , the two-body seagull amplitude S^(E 11 7 ) is given by a two-loop 
integral over the nucleon's momenta (see Fig. Wb): 
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<A',m'|S[ 2 ](£ 7 ,e 7 )|A,m> = jj ^<^( P ' - \*) *™ m2 (p - ~k 



(2 

x (m / 1 ,m / 2 |e-*e i S'] J 2] (-k / ,k;p / -k / ,-p';p-k, -p)\m 1 ,m 2 ). (5.3) 
To evaluate this amplitude, we perform integrations over p, p' and sum over spin variables 



numerically. Some details of the integration procedure are given in Ref. [46 




FIG. 7. The one- and two-body seagull amplitudes of 7c? scattering. 

In order to calculate the resonance amplitude _R(i? 7 ,0 7 ), we introduce the off-shell T- 
matrix of NN scattering, Tnn(E), and write the propagator G(E) = (E — H+iO)' 1 standing 
in Eq. ( |2.4| ) in the form 

G(E) = G {E) + G {E) T NN (E) G {E), (5.4) 

where Gq(E) = (E — H + iO) -1 is the propagator of free nucleons. Then R(E 1 , 9 7 ) turns 
out to be the sum of two terms, without and with NN rescattering in the intermediate state 
(see Fig. |8|; cf. Ref. ||38|| ). The term without rescattering reads 



(A',m'| J R™ t (E 7 ,6 7 )|A,m) = j (j£ - E k - *0 



(2tt) 3 W 

x (A', m'\T jd (-k'; p; p d + k)|mi, m 2 ) (mi, m 2 \T ld (k; p; p d + k)|A, m) 

+ (e <-> e'*, jfe <-> -A;'). (5.5) 

Here T ld {k; p, P) denotes the amplitude of deuteron photodisintegration without the final- 
state interaction (see Fig. ^) at the relative momentum p of the intermediate nucleons and 
the total momentum P. The energy Ek in Eq. ( |5.5|) is equal to 

E h = u + Pl_ A _(P* + k ) 2 , (5.6) 
k AM 4M ' V ; 

where A is the deuteron binding energy. The deuteron momentum p d is equal to — k in the 
CM frame, and it is unchanged when the crossing transformation k <->• —k! is applied. The 
procedure of a computation of T ld was the same as in Ref. j|(J , and we refer to this paper 
for further details and comments. Here we note only that the two-body contribution to T ld , 



2 f dp 1 1 

(m l ,m 2 \T\ 2 }{k;p;0)\X,m) = j — — *^ m , (p - -k 



(27T) 

x (mi,m 2 |e • j [2] (k;p, -p;p'-k, -p')\m' x ,m! 2 ), (5.7) 
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contains a loop integral over p'. Since we do not carry out angular integrations analytically, 
the evaluation of Eq. fl5.5|) actually involves a 9- dimensional numerical integration. Such 
a computational work was hard but it was done with reasonable computer resources. We 
have carefully controlled that the number of chosen nodes of integrations was sufficient to 
predict the observables of jd scattering like the differential cross section with the numerical 
accuracy better than 1%. 




a) b) 
FIG. 8. The resonance contribution R. Shown are terms without and with 7ViV-rescattering. 
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FIG. 9. Structure of the —* NN vertex. 
The resonance amplitude with NN rescattering has the form 

(A', m'\R—\E„ 6 7 )|A, m) = -J - E k - zO)" 1 -E k - zO)" 1 

x (A', m'\T jd (-k' ;p';p d + k)\m[,m' 2 ) (m[,m' 2 \T NN (E k ;p',p)\m 1 ,m2} 

x (m 1 ,m 2 \T jd (k; p; p d + k)|A,m) + (e <->• e'*, k <->• —A;'). (5.8) 

Here the A^iV scattering T-matrix is determined by the NN potential V through the 
Lippmann-Schwinger equation T NN (E) = V + VG (E)T NN (E). It is difficult to evalu- 
ate Eq. (|5.8|) straightforwardly. In order to simplify the computation, we used a separable 
approximation to T/vat. Actually, we took TV at from Ref. [|B8|j , in which the separable T- 
matrix was built for the Paris potential |89]]0 (see Ref. p5| for an explicit form of T^n given 
in our notation and normalization). 

Since the off-shell properties of Tjvtv for the Paris and Bonn potentials are not the same, 
the use of the "Paris" T-matrix in Eq. (|5.8| ) spoils the self-consistency of the whole calculation 
and even violates the gauge invariance of the resulting Compton scattering amplitude. A 
small mismatch appears between the resonance and seagull amplitudes at all energies, so 



To our knowledge, a separable approximation to the T-matrix of OBEPR does not exist in the 
literature. 
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that the balance ( 2.15|) prescribed by the low-energy theorem is not exactly fulfilled. At zero 



energy and forward scattering angle we get T s . a .(0, 0) = —0.47 for the spin-averaged (s.a.) 
amplitude of deuteron Compton scattering instead of the correct value of —0.50 (it is given 
in units of e 2 /M used here and below for the amplitude). For the spin-flip (s.f.) 'yd — > 'yd 
transition |1, —1) — > | — 1, 1) we get T s .f.(0,0) = —0.02 instead of zero. 

We believe, however, that the use of the "Paris" T-matrix in Eq. ( |5.8|) does not lead 
to practical problems at energies of a few tens MeV, because the rescattering contribution 
decreases with the energy. For example, at E 1 = 50 MeV and 7 = the amplitude (p~.8|) 
is equal to R r s e l cat = 0.009 — z0.076, what should be compared with the rest amplitude 
T s n ° rescat = —1.097 + z0.315 (found without nucleon-polarizability corrections). At larger 
angles the role of the rescattering contribution is even less important, as will be illustrated 
in the next Section. 

Given the Compton scattering amplitude {\' ,m'\T\\,m) , we find the differential cross 
section of jd scattering in the CM frame as 

7 A'm', Am 



where = y (2M — A) 2 + oo 2 and W = Ed + lo are the deuteron and total energies in the 
CM frame. The photon beam asymmetry is 

' ™'iTh ™\/v m 'iTi_i ™\* 

(5.10) 



da 1 - + daW 



2 Re E A 'm', m (^>' 


T 1, m) (A', m' 


T 


— 1, m}* 


J2x'm',Xm m ' 


T\\,m)\ 2 



where da L and dcr* are the differential cross sections for the incoming photons polarized 
perpendicular or parallel to the reaction plane xz, respectively. 



VI. RESULTS AND DISCUSSION 



A. Zero-range limit 



Before discussing results of the full model, let us consider the limiting case of a very weak 
binding of the deuteron. Specifically, let us assume that the typical nucleon momentum 
a = V MA, which is 45.7 MeV in reality, is much less than the pion mass. Said differently, 
we assume that the inter-nucleon distance r ~ a -1 in the deuteron is much larger than the 
AA-potential range. Moreover, we assume that the photon energy is also small, u = 0(a), 
so that all effects related with recoil corrections 1/M can be safely neglected too. In this limit 
the two-body contributions to the electromagnetic current and seagull become negligible, and 
so does the rescattering contribution ( p.8| ). Therefore, the Compton scattering amplitude is 
determined by the one-loop diagrams involving the operators S^, j' 1 ' and the asymptotic 
wave function of the deuteron, 



*ZrnM^^^^ m2 . (6.1) 
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Here C 1 md i is the Clebsch- Gordon coefficient. Keeping the terms of leading order over 

1/M in the electromagnetic operators and in the energy (|5.6|) of the intermediate nucleons, 
and calculating analytically the integrals ( |5.2| ) and ( |5.5| ), one arrives at the following 7c? 
scattering amplitude [j35| , [42ll :P| 

T(E J7 9 7 ) ~ ■ e'*{-F (g) - A + [(1 - «, - *0) 3 / 2 + (1 + W ) 3 / 2 ] }, (6.2) 



where 



4a q 

Foiq) = — arctan — (6.3) 
q 4a 

is the deuteron form factor in the considered zero-range limit, q — |k — k'|, and w = E y /A. 
The term with F (q) in Eq. ( |6.2|) represents the seagull contribution, whereas other pieces 
come from the resonance amplitude R(E y , 7 ). 

As a by-product of this computation, the total cross section a^T* vn {E 7 ) of deuteron pho- 
todisintegration can be derived through the imaginary part of the (spin-averaged) forward 
scattering amplitude ( |6.2|) : 



ImT s . a .(£ 7 , 0) = E,al d r pn (E,) * ^~ 2 (w - I) 3 / 2 . (6.4) 

When the photon energy becomes much higher than the deuteron binding energy A, 
the amplitude ( |6.2|) becomes equal to the proton Thomson scattering amplitude times the 
deuteron form factor F (q). This is just the seagull contribution S(E 7 , 6 7 ) to the amplitude 
(|6.2D. The rest (w-dependent) terms in Eq. (|6.2j) give the resonance amplitude i?(i? 7 ,6 7 ) 
which vanishes in the limit of E 7 ^> A. An instructive feature of Eq. fl6.2j) is however that 
this vanishing is rather slow, like oc E^ 1 ^ 2 . Therefore, the resonance amplitude can give a 
10—20% contribution to the differential cross section of 7c? scattering at E 7 ~ 100 MeV. 

In the opposite limit of very low energies, the binding corrections become large, and the 
amplitude ( |6.2|) is equal only one half of the proton Thomson scattering contribution: 

r(0,e 7 ) = -^e.e», (6.5) 



in exact agreement with the low-energy theorem for photon-nucleus scattering, Eq. (|2.15 ) 



The deviation of the nuclear amplitude (|6.5|) from the nucleon one, is related with the 



e 2 

resonance contribution -R(0, 7 ) which is equal to +^^ e ' e '* ^ n ^ ne considered zero-range 



12 Bethe and Peierls [35] who gave the very first analysis of jd scattering considered so low energies 
u ~ A <C a, at which the dipole El approximation is applicable. Accordingly, they had the form 
factor Fo(q) to be equal to 1. Equation (|6.2| ) as it is written here was given by Chen et al. [42]. 
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approximation. In the real case of the NN interaction of a finite range, both the seagull 
and resonance amplitudes get considerable modifications. For example, the resonance con- 
tribution i?M from the one-body electromagnetic current becomes smaller than e 2 /2M at 
zero energy ||90|| . Moreover, it depends on the deuteron spin. Omitting the rescattering 
correction (|5\8|) and evaluating Eq. (|5.5|), we have found 



e 2 

(±l,l|i2M noracat (0,0)|±l,l) = (±1,-I|i2 [1]norescat (0,0)|±l,-1) = 0.448—, (6.6) 



but 



p 2 

(±l,0|i? [1]norcscat (0,0)|±l,0) = 0.392— (6.7) 
for the Bonn potential OBEPR. 

e 2 

In contrast to w x \ the total resonance amplitude -R(0, B 7 ) is greater that ^j^ e ' e '* Du t 

it is also spin dependent. These features are easily seen from the relation ( |2. 15|) and from 
the explicit expressions ( p,17| ) and ( |2.18| ) for the one- and two-body seagull amplitudes at 
zero energy. 



B. Two-body seagull amplitude: low-energy behavior 

The two-body seagull contribution to the total Compton scattering amplitude dom- 
inates the binding corrections at energies of a few tens MeV, although the resonance contri- 
bution R is not negligible either. One can get more insight into a physical meaning of 
considering its low-energy behavior. Note that is a regular function of the photon energy 
below pion threshold, and it can be expanded in powers of uo. We have found that keeping 
terms up to order O(oo 2 ) is generally sufficient for getting quite an accurate approximation 
to results obtained through a numerical evaluation of Eq. ( p.3| ) at all energies up to 100 
MeV. The only exception is the contribution of the A-isobar, which requires also a 0(uj 3 ) 
term linear in the deuteron spin S. 

The spin-averaged part of the seagull amplitude at the considered energies can be 
described by the following most general form containing terms up to 0(u 2 ) and compatible 
with the discrete symmetries of the Compton scattering amplitude: 

+ AixA AoV<% + 4nA A/3w 2 (k • - kity). (6.8) 

Here, in our case, the number of neutrons, protons, and nucleons is equal to N = Z = 1 
and A = 2, respectively. The sign = is used in order to indicate that we have omitted pieces 
vanishing in the radiation gauge ( p,16|) . The four coefficients entering Eq. ( |6.8|) were found 
numerically to be equal to 
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« = 0.47, (r 2 K ) = 0.49 fm 2 



(6.9) 



and 

Aa = -0.72 x 10" 4 fm 3 , A/3 = 0.27 x 10" 4 fm 3 . (6.10) 

The coefficient k which determines the two-body seagull amplitude at zero energy is the same 
quantity which characterizes enhancement in the well-known electric-dipole photonuclear 
sum rule (by Thomas-Reihe-Kuhn, or TRK). Experimental data on the total cross section 
of deuteron photodisintegration seem to suggest somewhat lower value for k than that found 
from the Bonn potential, Eq. ( |6.9| ), namely k = 0.35 ± 0.10 (see, e.g., Ref. |HJ). For the 



OBEPR(B) version of the Bonn potential the difference would be even bigger since then n 
is predicted to be equal to ft = 0.50. 

Discussing a comparison with the experiment, it is worth to say that k as defined by Eq. 
( |B~8|) is not a direct observable but rather a (important) theoretical quantity which appears 
in the theoretical formalism after eliminating explicit meson degrees of freedom. Moreover, 
k may depend on the formalism specifically chosen, because the very separation of the total 
amplitude T into the resonance and seagull parts is subject to ambiguities. This, in turn, is 
because the separation of the Hamiltonian Q2.1|) into pieces of different order in the vector 
potential generally depends on the representation chosen and it is subject to unitary 
ambiguities (including off-shell ambiguities). In this respect, k in Eq. ( |6.8|) is as ill-defined 
and representation-dependent quantity as, for example, meson-baryon form factors which, 
nevertheless, are useful and meaningful theoretical objects. 

From a practical point of view, the representation ambiguity in the resonance and seagull 
amplitudes is probably not crucial. An often-used and convenient choice for fixing the 
resonance amplitude is to have R(u, 6) to vanish at high energies, as it is hold in the 
dispersion theory of nuclear Compton scattering ||. Yet, such vanishing is not exactly the 
case for our theory - in part, owing to the spin-orbit contribution and also because our 
theory is not relativistic. We refer to Ref. |J for further details concerning the relevance of 
the latter point. 

As for the above-quoted "experimental" estimate of k, it appears after two theoretical 
assumptions: (i) the validity of the so-called Gerasimov's argument |92| which makes a 
connection between the TRK and Gell-Mann-Goldberger-Thirring (GGT) dispersion sum 
rules having deal with the unretarded El and retarded total cross sections, respectively, and 
(ii) a little bit arbitrary cutoff at about 140 MeV in the GGT sum rule. Both assumptions 
are not strict. For instance, owing to a disbalance between corrections of higher order 
in v/c (viz., retardation and higher-mult ipole contributions), the Gerasimov's argument is 
actually violated in many models (see, e.g., Refs. p3| , |94| ,|6| for further detail and references). 



The quoted uncertainties in the "experimental" estimate for k represent only the cutoff 
dependence of the GGT dispersion integral, and one should take this into account when 
compares the "experimental" and theoretical predictions for k. 

The presence of the parameter (r 2 .) in Eq. ( |6.8| ) implies that the energy-independent part 
of the seagull has generally a g-dependent form factor, which is reduced to a linear function 
of the momentum transfer squared at low energies. Numerically, k is almost fully determined 
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by the pion exchange, Eq. ( [4.11] ), which gives k = 0.44, thus leaving only k hm = +0.03 for 
the contribution from heavier mesons of the Bonn potential. In contrast, the pure pion- 
exchange leads to a very small radius, (r 2 ) = 0.13 fm 2 . Actually, the most part of (r 2 ) 
comes from the pion exchange accompanied with the A-resonance excitation, as described 
by Eqs. ( ^20|) and flPip . 



The parameters Aa and Aj3 in Eq. (|6.8| ) determine the energy-dependent part of the 
seagull amplitude. Compared with Eq. (|3.6|), they can be loosely interpreted as medium mod- 
ifications to the electric and magnetic polarizabilities of the bound nucleon in the deuteron 
due to meson-exchange effects. Such quantities were introduced in this context by Hiitt and 
Milstein [[47]] (following the previous works 0,^]) and analyzed for spinless nuclei. See 
also Ref. [0], in which a review of a related experimental work is given. Similarly to k, the 
parameters Aa and Af3 are representation dependent and not direct observables. They are 
rather useful theoretical quantities which appear in the formalism with eliminated meson 
degrees of freedom. 

There is some distinction in the way how the quantities Aa and A(3 and the free-nucleon 
polarizabilities and /3n enter to the 7c? scattering amplitude. The medium modifications 
to the polarizabilities clearly have a non-local, i.e. two-body (and generally, many-body) 
origin, and they are expected to be accompanied with a "two-body" form factor F 2 (q) which 
describes a distribution of the center of the relevant nucleon pairs in the nucleus (in the 
deuteron we would expect F 2 (q) = 1). The form factor F 2 (q) should be different from 
the usual "one-body" form factor F(q) describing the distribution of single nucleons in the 
nucleus and accompanying the contributions of the free-nucleon polarizabilities. 

In the deuteron case the difference between F 2 (q) and F(q) is especially large since the 
radius squared of the one-body form factor (averaged over the deuteron spin) is quite large: 
(r 2 ) = 3.9 fm 2 . See Refs. [87, fj for a more quantitative analysis of this difference in the 



case of heavy nuclei.0 Numerically, the values ( |6.10| ) are dominated by the (retarded) 
pion exchange which gives alone Aa = -0.99 and A/3 = 0.36 (in units of 10" 4 fm 3 ). The 
retardation effects incorporated through Eq. ( |4.23|) are very important here, and they give 
alone Aa( ret ) = -0.84. The values (|5TTUD are similar but essentially larger, especially for 



Aa, than estimates obtained in Refs. |$7,|] for the lightest even-even nucleus 4 He on the 
basis of the correlated Fermi-gas approximation which is suitable for heavy nuclei. 

Considering the spin-dependent part of the seagull amplitude in the similar way, one 
has to introduce a few more parameters. We will not discuss all of them here and mention 
only two, the tensor enhancement parameter kt and the tensor modification of the electric 



13 The presence of the radius (r 2 ) in the two-body contribution to the seagull amplitude ( |6.8| ) 
implies that there is no universal two-body form factor which accompanies both the energy- 
independent and energy-dependent part of the seagull. Due to the term with (r 2 ) > 0, the 
form factor which multiplies the energy-independent part (i.e. k) has a larger radius than that 
multiplying the energy-dependent part. This feature was also first found in Ref. [87] for heavy 
nuclei. 
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polarizability Acty. They appear in the low-energy expansion of the tensor part of S\j :0 

X 



SiSl + SlSi -\ S(S+mi ] + ... (6.11) 



We have found numerically that k t = 0.24, so that the seagull amplitude at low energies 
has a strong spin dependence. This number is again dominated by the pion exchange. The 
heavier mesons of the Bonn potential OBEPR give only k^ m = —0.03. As for Aa^, it 
gets the largest contribution from the retardation effects in the exchange-pion propagator 
which give alone Aa { * et) = -1.3 x 10 4 fm 3 . Therefore, in contrast to the case of heavy 



nuclei considered by Hiitt and Milstein [17|, the meson-exchange-induced modification of 



the electric polarizability of the bound nucleon is essentially deuteron-spin dependent. 

The retardation correction which manifests itself in the parameters Aa and Aolt in- 
creases noticeably the differential cross section. For example, this increase is equal to 5—7% 
at 100 MeV at all scattering angles. 

The energy dependence of the two-body seagull contribution at forward scattering 
angle is illustrated in upper panels of Fig. [10]. It is rather flat in the case of the spin-averaged 
amplitude, 

SSL{E y ,0) = \ £ (l,m|^(£ 7 ,0)|l,m), S£L(0,0) = -^*, (6.12) 

■J m=-l,0,l zm 

and it is more pronounced in the case of the spin-flip amplitude, 

53.(E 7 ,0) = (-1,1|^(E 7 ,0)|1,-1), S [ 3X0,0) = ^k t . (6.13) 

The A-resonance contribution becomes rather noticeable above 60 MeV, and it diminishes 
the energy dependence introduced by the pion exchange. As the result, the energy de- 
pendence of the seagull amplitude S^(E 7 ,0) gives only a 2% increase (in the absence of 
the polarizability contribution from the free nucleon) in the differential cross section of 7<i 
scattering at forward angle and E 1 = 100 MeV. 

In the bottom panels of Fig. |H] we show the spin- averaged and spin-flip two-body seagull 
amplitudes in the case of backward scattering. They are defined by equations like (|6.12|) 
and QS.13| ), in which 7 = is replaced by 7 = 180° and the helicity of the final photon 
is inverted. We see that the energy dependence of the spin-flip amplitude is rather steep in 
this case. It increases the differential cross of deuteron Compton scattering at 100 MeV by 
5% (again, in the absence of the free- nucleon polarizabilities). 



14 A complete basis for representing spin-dependent Compton scattering amplitudes in the general 
case of spin S > 1 was found by Pais pi . 
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FIG. 10. The energy dependence of the spin-averaged two-body seagull amplitude SLL and 
the spin-flip amplitude S^j . Units are e 2 /M. The (retarded) pion-exchange contribution is shown 
in dotted lines. Successive additions of heavy mesons and the A-isobar leads to the dashed and 
solid lines, respectively. 



In general, we have found that the effect of the A excitation onto the seagull amplitude 
S^ n is not large. It does not exceed —2% in the differential cross section of 7c? scattering, 
da/dfl, at all considered energies. The two-body contribution of the A-isobar to the reso- 
nance amplitude R was found to be not large too. For example, it changes the differential 
cross section at 100 MeV by +3% at forward angle and by —0.6% at backward angle. At 
this point we disagree with the results of Ref. |37 |, in which it was found that the effect of 
the A-excitation onto the resonance amplitude at 100 MeV is negligible at small angles and 
very large at 7 > 90° giving a +15% increase in da/dVL. 



C. Model dependence 



Contributions of the different components _R norescat ; gPl and R rescat of the total 
amplitude T to the differential cross section at a few selected energies are shown in Fig. |Tl~l.Pl 
It is seen that the effects of the resonance amplitude, the one-body seagull as well as the 



15 Since the seagull and resonance contributions are not gauge-invariant separately, we remind that 
we use the gauge ( 2.16|) to calculate R and S. 
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two-body seagull are of similar scale, though the total effect of the seagulls is more than 
70% in the energy region under consideration. At the same time rescattering has a modest 
impact on the differential cross section. Our results confirm findings of other approaches 
p8|^40| that the rescattering decreases the forward differential cross section and that this 
decrease is between 7% to 12% in the energy region of 50 to 100 MeV. 

There was some discrepancy in the previous work concerning the role of R Tescat at back- 
ward angles. Weyrauch |3S| found that the rescattering does not contribute at backward 
angles at all, whereas in the later calculations |39|,f|(J a visible increase in da/dQ was claimed. 
Our results agree with the latter conclusions and suggest that the effect of i? rescat (i? 7 , 180°) 
ranges between +7% at 50 MeV to +3% at 100 MeV. Since an accurate calculation of _R rescat 
is a difficult problem for our computational scheme, this finding of a relatively small effect 
of _R rescat at energies and angles where experimental data are available provides some justifi- 
cation to our approximate use of a separable potential |88| for a computation of the off-shell 
T-matrix of NN rescattering. 
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FIG. 11. Contributions of different parts of the amplitude to the differential cross section (CM) 
of jd scattering at 50, 70, and 100 MeV. The contribution of the resonance amplitude Jl norescat 
alone is shown in dotted lines. Successive additions of the one-body seagull the two-body 
seagull S™, and the rescattering amplitude _R rescat give the dashed, dash-dotted and solid lines, 
respectively. Nucleon polarizabilities are not included into 5W. Data are from Ref. [99]. 



An instructive feature of the calculation is that the spin-orbit (s.o.) contributions to the 
electromagnetic current and seag ull, j [1]s o ' and SW s -°-, which are relativistic corrections of 
order 1/M 2 in Eqs. (|3.5|) and ( |3.6|) , are rather essential. Previously, it was found that 
the spin-orbit current is responsible for the long-standing discrepancy between the theory 
and data on deuteron photodisintegration at forward and backward angles. In the reaction 
of Compton scattering, the spin-orbit interaction is important not only at extreme angles, 



and this importance increases with the photon energy. Either j 



[1]b.o. 



or 



S [1]s -°- leads to an 



approximately equal decrease in the differential cross section at the forward angle, and the 
total effect of the spin-orbit interaction on da(E 1 ,0)/dQ is —4% at 50 MeV and —15% at 
100 MeV. A somewhat different situation happens at the backward angle. The spin-orbit 
current jW s '°- still decreases the differential cross section, but the spin-orbit seagull S'W 8 ' ' 
makes a bigger increase. The net effect of the spin-orbit interaction at 180° is +4% at 
50 MeV and +8% at 100 MeV. In the central angular region of 6 7 ~ 90°, the spin-orbit 
interaction has a little impact on da/dQ ranging between +0.8% at 50 MeV and —0.2% at 



33 



100 MeV. 



Staying within the Bonn-potential picture, we have checked how the differential cross sec- 
tion depends on a specific choice of the potential's parameters. The OBEPR and OBEPR(A) 
versions of the Bonn potential give da/dfl which are different at most by 1% in the energy 
range of 50—100 MeV. A bigger difference is found for the OBEPR and OBEPR(B) ver- 
sions, though it decreases with the photon energy. For example, the OBEPR(B) potential 
gives da/dVL which is bigger by 5% (7%) at 50 MeV and 0.5% (5%) at 100 MeV for for- 
ward (backward) angles, respectively. The main reason for such a difference comes from a 
very large value for the cut-off parameter A„- = 2 GeV used in the 7r-exchange potential 
of the OBEPR(B) version, whereas A„- = 1.3 GeV for OBEPR. Respectively, the seagull's 
parameters k and are bigger for OBEPR(B) too.Q 

We may note that the OBEPR(B) version of the Bonn potential does not provide a 
satisfactory description of observables in deuteron photodisintegration jKJ and thus it is 
not fully realistic. Therefore, one may conjecture that the sensitivity of the results on '-yd 
scattering would not be so noticeable if one restricts oneself to "realistic" potentials only. In 
this respect it is worth mentioning that the use of different momentum-space versions of the 



Bonn OBE potentials was found to yield version-independent results for <yd scattering 
within 1%. 



The present results are in a qualitative agreement with our previous calculation |40fl done 
in the framework of a "minimal model", in which MECs and MESs are evaluated through 
the minimal substitution p — > p — eA in the (Paris) NN potential. As they were published, 
those older results did not include the spin-orbit interaction. After taking into account jW s -°- 
and 5 , W S "°' the predictions of the minimal model become closer to the results of the present 
work, especially as for the shape of the angular dependence. Still, the minimal model gives a 
lower differential cross section: by 6% at 70 MeV and by 11% at 100 MeV. Such a difference 
can be traced in part to a weaker energy dependence of the seagull amplitude found in the 
minimal model and to the absence of the two-body A-isobar effects in that model. 

A comparison of the present predictions with the results of other calculations of 80 's - mid 



90 's |37|- |39fl is shown in Fig. [12|. There is a reasonable agreement between all the predictions 
at low energies like 50 MeV, maybe with the except for those of Ref. [|37[| . When the energy 
increases up to 100 MeV, we predict a bigger angular variation of da/dfl that other works 
do. In fact, the discrepancy between all the results of different authors is dramatically large 
at 100 MeV. Moreover, since the spin-orbit interaction was not taken into account in Refs. 



37H39| and since the spin-orbit effects decrease the forward cross section by 15% at 100 



MeV, the genuine disagreement of other components of the scattering amplitude with those 



of Ref. [39] is even more serious than Fig. 12 suggests. 



16 As a word of precaution we have to remind that we always calculate the rescattering contri- 
bution ij rescat using a fixed T-matrix obtained with the (separable-approximated) Paris potential. 
Therefore, the specific numbers indicated in the above discussion do not show the full change 
in the theoretical predictions when the potential is changed, although we assume that they are 
qualitatively correct. 
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FIG. 12. Differential cross sections (CM) of 'yd scattering. Dotted lines: Weyrauch and 



Arenhdvel [37]. Dashed lines: Weyrauch f38|| . Dash-dotted lines: Wilbois et al. [39|. Solid lines: 



the present work. Nucleon polarizabilities are turned off. Data are from Ref. [99]. 



It is interesting that the very low differential cross section obtained by Wilbois et al. 
at 100 MeV has found a support from a recent work of Karakowski and Miller (KM) 



see Fig. [13|. Below we mention a couple of possible reasons why the KM model disagrees 
with our model and shows a wrong behavior at high energies E 1 ~ 100 MeV.Q First, 
we have rather different effects of the spin-orbit interaction. In Ref. HJL this interaction 



was taken into account partly, i.e. only through the electromagnetic seagull, not through 
the electromagnetic current. However, it gave a much bigger decrease in the differential 
cross section at forward angles than that we found (cf. Fig. 12 in Ref. ||41|| ). The second 
reason, perhaps less important numerically, might be that there is a mismatch between 
the electromagnetic and strong-interaction parts of the Hamiltonian i^KM of the KM model 
which destroys the gauge invariance and actually signals that some electromagnetic charges 
or currents in the system are missing in the theoretical formalism. The mentioned mismatch 
is that the electromagnetic two-body part of Hkm includes only the point-like pion-exchange 
piece, whereas the wave function of the deuteron is constructed using a more complicated 
(and more realistic) Bonn potential.^ The violation of the gauge invariance in the KM 
model did not lead to visible problems at low photon momenta fcr < 1 owing to the use 
of the Siegert transformation which ensured automatically the fulfillment of the low-energy 
theorem Q2.15 ). However, when kr becomes large (this is the case for energies E 1 ~ 100 



MeV), the Siegert transformation does not help, and the missing charges and/or currents 
can become important. 

There is a simple way to verify the theoretical calculations in the particular case of 
7 = and to see that the predictions of Ref. pO at high energies are invalid. Using the 



We are indebted to Prof. G.A. Miller for useful comments made on this point. 

18 Actually, the procedure of Ref. [O] is even more complicated in this respect, because one more 
potential (Reid93) is used in another part of their computation, when the rescattering correction is 
found. To be honest, we have to say that we also do a similar mixture in order to evaluate _R rescat 
(see Section |V|) . 
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FIG. 13. Differential cross sections (CM) of scattering. Dashed lines: Karakowski and 
Miller Dotted lines: Beane et al. @. Dash-dotted lines: Chen et al. |||]. Solid lines: the 
present work. Nucleon polarizabilities are included, and ajy — (3n = 9 is used for drawing the solid 
curves. Data are from Ref. p9] (solid circles) and Ref. |102(| (open circles). 



Gell-Mann-Goldberger-Thirring dispersion relation jointly with the optical theorem (cf. Eq. 
(|6.4|)) for the spin- averaged amplitude T s . a ., we write 



ReT s . a .(w,0) 



AM 7T Jo ,c 



,2 



doj' 



or 



(6.14) 



where <7 to t is the total photoabsorption cross section. Keeping in mind that the total cross 
section of meson production off the deuteron is dominated by meson production off quasi- 
free nucleons, we see that this part of the photoabsorption cross section is responsible for the 
component of the scattering amplitude which is related with the polarizabilities of free 
nucleons (up to relatively small effects due to medium modifications of these polarizabilities). 
Therefore, subtracting the meson-production part of the photoabsorption cross section and 
keeping in Eq. ( |6.14| ) non-mesonic, or photodisintegration part of the cross section, we can 
approximately identify the resulting r.h.s. of Eq. ( |6.14| ) with the jd scattering amplitude, in 
which the internal (mesonic) structure of the nucleon is disregarded. In other words, using 
the deuteron photodisintegration cross section cr^t^ pn instead of cr tot in Eq. ( |6.14|) , we should 
obtain the 7c? scattering amplitude with point-like nucleons having zero polarizability. See 
Ref. for a more detailed discussion of these steps. 



We have evaluated the integral in Eq. ( |6.14| ) at energies u < 100 MeV using: (i) the 



effective-range parameterization of <Jt£~* pn (u') at energies uj' below 20 MeV (see Ref. ||100 
Eq. (2.18)) which gives an accurate description of experimental data at low energies, and 
(ii) a phenomenological fit ||101|| to available experimental data between 20 and 440 MeV. 
At higher a/, the photodisintegration cross section is small and can be safely neglected in 
the integral. The result of such an evaluation of Eq. ( |6.14|) is shown in Fig. ([[4]) by the solid 
curve, together with our predictions (dashed line) based on the Bonn-potential picture, in 
which the nucleon polarizabilities are disregarded. Generally, we find very good agreement 
between the two curves. Some disagreement of about 6% at very low energies appears due to 
an approximate way of finding the rescattering amplitude _R rescat ; as it was already mentioned 
111 Section 0. At energies above 10 MeV, the rescattering amplitude is less important, and 
the agreement between the two calculations improves. It is better that 3% even at 100 MeV. 
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It is needless to say that the Bonn-potential picture nicely reproduces the experimental data 

at all energies below pion threshold as well as the differential 



on the total cross section aj^ pn 



cross section of deuteron photodisintegration and polarization observables [[K 
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FIG. 14. The spin-averaged amplitude of forward 7d-scattering. Solid line: the r.h.s. of Eq. 
( 6.14| ). Dashed line: results of the present work based on the Bonn potential. Meson photopro- 
duction and nucleon polarizabilities are disregarded. 

Now let us compare our predictions with those obtained within two different flavors of 
effective field theory (EFT) for few-nucleon systems @,fl4j] ( a general review of the EFT 
approach to nuclear problems can be found in Ref. [ |103|| ). At "high" energies E 1 ~ 70—100 
MeV, we have a qualitative agreement with the results of Beane et al. |44], who used the 
so-called Weinberg formulation of the nuclear EFT. See Fig. [l3| for a comparison. In part, 
a proximity of our and Beane et al. predictions is caused by their use of a realistic wave 
function to evaluate matrix elements of the ChPT kernel, the latter denoting the amplitude 
of syNN — > syNN taken to order 0(Q 3 ) in the chiral perturbation theory expansion. It 
was actually the wave function of the Bonn potential.^ Furthermore, the dominating part 
of the two-body seagull operator is the same in our two approaches. It comes from the 
retarded pion exchange. Despite the Bonn-potential picture uses form factors in the 7rNN 
vertices and additional heavier mesons (to make improvements to the potential at small 
distances), the effects of the form factor and the heavy mesons onto the seagull amplitude 
are reduced to some changes of k which almost cancel each other. Specifically, we have 
found ft = 0.48 — 0.04 + 0.03, where the three numbers are the contribution of the pure pion 
exchange, the contribution from the fictitious A particle simulating the irNN form factor 
(see Section |TV|), and the contribution of p, u>, 5 and a mesons, respectively. 

H 



We can notice that all the contributions counted in Ref. 



have been taken into account 



in our calculation as well. Beyond that, we included other corrections which formally belong 



19 The optional use of the Nijmegen-potential wave function (see Fig. 15 in Ref. ^4|) leads to 
bigger deviations from the dispersion predictions and from the Urbana data |99|| . It was claimed 
in Ref. @ that any wave function of the deuteron with the correct binding energy, thus including 
for instance ^(p) given by Eq. (6J), could also be used within the considered order 0(Q 3 ) of the 
power expansion. Being valid theoretically, this makes, however, a big numerical difference! 
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to higher orders in Q in the power counting scheme of the EFT but are large numerically. 
An instructive example is the amplitude R norescat given by Eq. ( |5.5|) , in which each of the 
amplitudes of the 7c? — > pn transitions is dominated, at high energies, by the pion-exchange 
current j^ 77 , as shown in Fig. [| This rather large contribution formally appears only in 



order 0(Q 5 ) of the power counting scheme of Ref. [pB |. 

A resemblance of the two predictions seems to be lost at lower energies, where the 
scattering amplitude of Beane et al. begins to deviate from the correct value fixed, for 
instance, by the dispersion relation (|6.14|) at the forward angle. Such a failure is not a surpise 
and it was anticipated in Ref. as a result of contributions from the NN intermediate 
states which break down the used power counting at low energies. 

It is worth mentioning that the evaluation of matrix elements of the ChPT kernel between 
phenomenological (Bonn or Nijmegen) wave functions, which are not consistent with the 
one-pion-exchange dynamics of the NN interaction incorporated into the ChPT kernel, 
automatically means the absence of the gauge invariance in the scheme of Ref. MM. This 



itself is a sufficient reason for a failure of such a theory at very low energies where the gauge 
invariance is crucial. 



we can 



In view of close magnitudes of k arising in our approach and in that of Ref. [4"3 
conjecture that the main difference between the two predictions at energies of about 100 
MeV is related with our taking into account the A-resonance excitation (both in one-body 
and two-body operators) and with our taking into account the contribution i^ norescat ( Q f the 
correct magnitude) and the contribution _R rescat . 

The problem with the gauge invariance and with the region of very low energies does not 
exist in the version of EFT used by Chen et al. [42]. Their work is based on the so-called 



Kaplan-Savage-Wise (KSW) regularization which successfully resolves the problem of a 
poor power-series convergence in the case of large s-wave AA-scattering lengths ||104jl . NN 
rescattering contributions are accurately taken into account in that approach. However, even 
being quite accurate at low energies, this approach becomes inapplicable when the momenta 
of nucleons in the rescattering diagrams exceed the range of convergence of power series 



which is about A^at = ^^M/ g^ NN ~ 300 MeV. This makes predictions of Ref. |42] not well 



controlled at energies > 70—90 MeV. Therefore, there is no surprise that these predictions 
at 70 MeV lie visibly lower than our predictions (and those of Beane et al. [PJ), including 
the 7 = point, where the dispersion relation ( |6.14|) strongly favors our calculation. It is 
worth mentioning that neither A-isobar excitation nor the spin-orbit current and seagull are 



taken into account in Ref. [42[ since these pieces appear only in higher orders of the used 
expansion. 



D. Determination of the nucleon polarizabilities 



Considering the nucleon dipole polarizabilities «jv and (5n in the electromagnetic seagull 
operator (|3.6|) as free parameters, we can check the sensitivity of the differential cross section 
da/dfl with respect to variation of these parameters. Our results are shown in Fig. |T5| 
together with a few experimental data available from Urbana (E 1 = 49 and 69 MeV) |M| and 
Saskatoon (E^ = 94 MeV) ||102|| . We do not show how the differential cross section depends 
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on the sum of the electric and magnetic polarizability, because this sum is reasonably-well 
fixed by the Baldin sum rule, Eq. ( p. 12 ). As for the difference of and (3^ which is not 
well-known theoretically, it can be determined from data at large scattering angle. 
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FIG. 15. Dependence of the differential cross section (CM) of yd scattering on the nu- 
cleon-averaged dipole polarizabilities a at and (3n- Dashed lines: all the polarizabilities (including 
those of higher order) are turned off. Dotted lines: only higher-order polarizabilities, Eq. fl3,13|) , 
are included. Dashed-double-dotted, solid, and dashed-dotted lines: a n — Pn = 6, 9 and 12, re- 
spectively, oat + (3n = 14.6 is fixed, and the higher-order polarizabilities are included. Data are 
from Ref. p9| (solid circles) and Ref. ||102j] (open circles). 



Of course, the highest sensitivity is observed at the highest photon energy. Nevertheless, 
we believe that data at medium energies like 70 MeV are also quite useful, because theoretical 
uncertainties in our computation related, for example, with omitted relativistic corrections 
or with omitted dispersion effects due to two-pion exchanges are expected to be smaller at 
lower energies. 

Dotted lines in Fig. 15 illustrate the fact that the higher-order polarizabilities (|3.13|) are in 
no way negligible when the nucleon dipole polarizabilities are determined from yd scattering 
at energies ^ 70 MeV. This feature was paid attention to also in Ref. [[Hj], in which the 
higher-order contributions appeared as an intrinsic part of the one-pion-loop diagrams of 
the ChPT kernel. Since, however, we use the higher-order polarizabilities which are given 
by dispersion relations P5]J52]Jr]^|J57[] and which are very different from those suggested by 
the one-pion-loop mechanism (in part, due to the A-contribution, see Ref. [p5|] ), we predict 
a much bigger effect at backward angles. 



A straightforward two-parameter fit of the Urbana data gives = 14.5 ± 2.7 and 
Bn = 6.6 ± 2.7, whereas a similar fit of the Saskatoon data | |102| | gives a lower value of the 
electric polarizability: = 8.4 ± 1.8 and (3n = 6.2 ± 1.8. Making a combined fit of all the 
data, we obtain 



a N + (3 N = 17.1 ±1.6, 



what is in agreement with the theoretical expectation ( |3.12| ), and 

~ Pn = 4.0 ± 1.5, 



(6.15) 



(6.16) 



though with a poor xV-^d.o.f. = 21/9. Systematic uncertainties of the experimental data 



are included into the obtained estimates (|6.15|) and Q5.16j ). However, it is not so easy to 
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estimate uncertainties introduced by the theoretical model. Certainly, they are not less than 
the experimental uncertainties. 

Taken as they are, these numbers, together with the experimental data on the polariz- 
abilities of the proton ( |1.2|) can be considered as an indication that the electric polarizability 
of the neutron is a n = 9 ± 3, and the neutron magnetic polarizability is (3 n — 11 ± 3. While 
the obtained sum a n + /3 n = 20 ± 3 reasonably agrees with the theoretical estimate (|3.11b|) , 
the obtained difference a n — (3 n = —2 ± 3 is rather far from both the similar difference in 
the proton case found experimentally, a p — (3 P = 10 ± 2 |L2 , 16| , and from theoretical esti- 
mates based on dispersion relations which predict roughly a n — j3 n ~ a p — /3 P (see, e.g., Refs. 
|I|j3|j6^] ) . It is clear that a further experimental and theoretical work is needed to reduce 
the uncertainties. New data can appear from Lund [ |105| ] . 

Among other observables of 'yd scattering which are sensitive to the nucleon polariz- 
abilities too, we briefly discuss the beam asymmetry E, Eq. ( [5.101 ). In Fig. |Tj] we show 
how different components of the Compton scattering amplitude affect E (this is helpful for 
imagining a possible scale of model uncertainties) and how E is sensitive to the nucleon 
polarizabilities. One can notice a strong dominance of the one-body seagull amplitude, 
whereas the role of the two-body seagull contribution is smaller than that in the case of the 
differential cross section. The role of NN rescattering is again small. 

The spin-orbit interaction essentially affects E and gives a 10% increase at central angles. 
It mainly comes through the one-body seagull amplitude S^. The contribution of the A 
excitation to the two-body seagull amplitude changes the beam asymmetry by less than 
1%, but the A contribution into the resonance amplitude R is rather visible, reducing E 
by 6% at 100 MeV. Pion- retardation effects have only a tiny impact on E. The use of the 
potential OBEPR(B) instead of OBEPR has a big effect and reduces E by 14% at 100 MeV. 

It looks like experiments with the linearly-polarized photon beam can also be useful for 
measuring the nucleon polarizabilities, provided the accuracy of measurements is better than 
~ 5-10%. 
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FIG. 16. Left panel: Contributions of different parts of the Compton scattering amplitude to 



the beam asymmetry at 100 MeV. Notation of curves is as in Fig. |11|. Right panel: Dependence 
of the beam asymmetry on the nucleon-averaged dipole polarizabilities a at and (3n- Notation of 



curves is as in Fig. 15 



We conclude saying that the reaction of deuteron Compton scattering at energies of about 
50-100 MeV has a great potential for a determination of the electromagnetic polarizabilities 
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of the neutron. Currently, the available theoretical models show a big divergence in their 
results, in part because they do not take into account all important contributions. So, a 
further theoretical work is needed to improve the accuracy of the models before any firm 
conclusions could be inferred about the values of ajy and (3^. Better experimental data are 
also needed to this aim. 
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APPENDIX A: ELECTROMAGNETIC SEAGULLS FROM HEAVY MESONS OF 

THE OBE POTENTIAL 



In this appendix we give explicit formulas for the seagulls Sfj produced by the meson 
exchanges a = 77, a, 5, u and p of the Bonn potential (OBEPR). All of them are obtained 
through a direct evaluation of the diagrams shown in Fig. |l]c. The electromagnetic effective 
meson-nucleon vertices jaNN and r yyaNN in these diagrams arise from the relativistic 
boson- nucleon effective Lagrangian of Refs. [|73| , [73f| , in which a nonrelativistic reduction is 
done and terms up to order 0(M~ 2 ) are only retained. Such a procedure is consistent 
with the whole construction of the OBEPR, because this potential itself is built through 
the truncation of the relativistic Feynman diagrams of the one-boson exchanges to order 
0(M~ 2 ) (see Ref. [Tj, Appendix A. 3). 

Technically, the nonrelativistic reduction can be conveniently performed [|8T| by consid- 
ering appropriate relativistic Feynman diagrams (see Fig. |I7p and keeping only the negative- 
energy part P- of the nucleon propagators. To leading order in 1/M, it is sufficient to take 
P_ in the static limit, i.e. 

Ml ?)■ 



It is worth noticing that the contact vertex jaNN in Fig. 17 appears only in the case of 



a = p, being caused by the tensor coupling of the charged p-meson to the nucleon. In 



the formalism of the pseudo-scalar irNN coupling used in Refs. |73| . |74]j , the contact vertex 
^/7tNN is absent. 



111 I II 

III I II 

a 1 1 1 a 1 11 

a) b) 

FIG. 17. Effective contact meson-nucleon vertices 'jaNN (diagrams a) and 'j^aNN (diagrams 
b) arising from the antinucleon degrees of freedom. Thick lines denote the negative-energy part P_ 
of the nucleon propagator, Eq. (|Al|). 

We use the same notation as in Section [TV]. In particular, the momenta q x , q 2 , Ki and 
K 2 are defined by Eqs. ( [4.8|) and (|4.12| ). We introduce also the vectors 



Pi = Pl + p[, P 2 = p 2 + p 2 . (A2) 

The functions G a (q), G\ a (q 1 , q 2 ), G , 2 Cf (q 1 , q 2 , K), and Z) Q ,(q 1 , q 2 , Ki, K 2 ) used below are 
defined as in Section [TV] with the evident replacement tc — > a in all masses m a and cut- 
off parameters A a . The constants n a denote the ratios of the tensor and vector coupling 
constants for the uNN and pNN vertices, n a = f a /g a . 
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Evaluating the diagrams in Figs. |T]c and [17|, we obtain the following results. 
Isoscalar exchanges: 

£*£jSij(- k '> k ; P'n P25 Pi> P2) = °. 



ef e^(-k', k, p[, p' 2 ; Pl , p 2 ) = -e • e'* |^| [^(q,) + (1 <- 2) 



e i e j^.(-k', k > Pi' P2! Pu P2) 



e ■ e 



2M 2 



^C^q,) - 2Z 1 Z 2 G U1 (K 1 ) + (1 <- 2) 



Writing the last equation, we have used that k u = for the Bonn potential. 
• Isovector exchanges: 



e '* e j S ij(- k 'i k ! Pl> P2; Pl> P2) = -c ■ 6 



2M 2 



ZnZGsiqJ + (1 <-> 2) 



+ Ia§{ ( q2 • ^"(^ Kl ) [ T ^(qi ■ e + ^1 x Pi • 6) 



+ i(ri x t 2 ) 2 (Pi • e + z<Ti x q x • e)J + (1 «-> 2) 
+ (e^e'*, Kx^-K 2 )} 

P 2 + P 2 -q 2 -q 2 



-e g 5 Ti 2J D 5 (qi, q 2 , Ki, K 2 ) 



16M 2 



.<Ti x q : • Pi + cr 2 x q 2 ■ P 2 



8M 2 



eV 



ef e^-k', k; P ;, p 2 ; Pl , p 2 ) = -e ■ e'* ^jj§ (ZiT 2 2 - Kp T 12 )G p (q 2 ) + (1 <-> 2) 



2„2 



+ 4^{G>( K i) [e ■ e'*m 2 - 4z\Z 2 ) + (1 + k p ) 2 T 12 CTl x e ■ <x 2 x e' 



+ 



4M 2 



- (1 + k„)(ti x r 2 ) z (<r a + <r 2 ) ■ e'* xe]+(ln 2)} 
(q 2 • e'*G lp (q 2) K!){T 12 [(1 + 4K p ) qi • e + + 2/c p )«n x Pi ■ e 
- 2i(l + Oo-i x P 2 ■ e - 2(1 + k p ) 2 (Ti x e ■ <r 2 x q 2 



+ z(ri x t 2 ) 2 (Pi - 2P 2 ) • e + z(l + 2k p )<ti xq r e- 2i(l + k p )ct 2 x q 2 
+ (1 <-> 2)) + (e ^ e'*, Ki <-> -K 2 ) 

P 2 + P 2 - 4Pi • P 2 



+ e g 5 T 13 D p (q 1 , q 2 , Ki, K 2 



1 + 



+ i(l + 2K f 



16M 2 

CTx x q 1 ■ Pi + cr 2 x q 2 ■ P 2 

8M 2 



(1 - 4k, 



qi + qj 
16M 2 
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ill + k c 



<j\ x q 1 ■ P 2 + <x 2 x q 2 • P 2 



+ (! + «, 



4M 2 

2 Q-! x qi ■ cr 2 x q 2 

j AM 2 



(A7) 



Writing Eqs. (|A6|) and flA7j), we used the radiation gauge (|2.16|) . 

With the help of Eqs. ( [4.9|) and ( |4.16| ), one can verify that thus constructed seagull 
operators Sfj satisfy the equation ( |4.17| ), provided the electromagnetic currents j a are taken 
as obtained [^6j from the same boson exchanges to order 0(M~ 2 ). 

Numerical values of the masses m a , the couplings g a (as well as n a for the vector mesons), 
and the cutoff parameters A Q for different bosons are taken exactly the same as for the Bonn 
potential (OBEPR) itself (7|tt The only exception concerns the cx-exchange. The matter 
is that the Bonn parameterization of the a-exchange suggests to use a different mass m a 
and the coupling g a for different NN channels with the total isospin 7 = or 7 = 1. We 
found this feature inconvenient for building electromagnetic operators which mix the isospin. 
Since we noticed no visible distinction between our predictions using <r-MEC and cr-MES 
with either of the two sets of (m CT , g a ), we took for the sake of simplicity the cx-meson 
parameters proposed by the Bonn group for the 1 = channel. 

As a final comment we have to mention that, strictly speaking, the nonrelativistic reduc- 
tion of the Feynman diagrams to order 0(M~ 2 ) considered as a method of a determination 
of the OBE potentials V a , the OBE electromagnetic currents j", and the OBE electromag- 
netic seagulls may need a further refinement. The matter is that the operators obtained 
in this way are manifestly frame- dependent. Specifically, they explicitly depend on the in- 
dividual momenta of the nucleons, Pi and P 2 , rather than on the relative variable Pi — P 2 
(see, e.g., Eqs. (|A6|) , ([AT]) and also Ref. [@). This is not what is expected for potentials, 
MEC, and MES in the nonrelativistic framework. 

It is possible to propose a modification of the diagrammatic representation of the boson 
exchanges which leads to Galilei-invariant results and to some changes in the above equations 



for MEC and MES ||106|1 . We checked, however, that such a modification has only a minor 
numerical effect and can be neglected in the present context. 
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